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(54) SYSTEM FOR EVALUATING PRICE RISK OF RNANCIAL PRODUCT OR ITS FINANCIAL 
DERIVATIVE, DEAUNG SYSTEM, AND RECORDED MEDIUM 



(57) A system for correct evaluating a price distri- 
bution arYd a risk distnl^ution for a financial product or its 
derivatives is provided by Introdudng a probability den- 
sity function generated with a Boltzmann model at a 
higher accuracy than the Gaussian distribution for a 
probability density. A dealing system applying this price 
evaluatk>n technique is also provided. 

The price and risk evaluation system has an initial 
value setup unit 3 and an evaluation condition setup unit 
4. Initial values include at least one of the price, the price 
change rate, and the price change directk>n of a financial 
product. The evaluation conditk>ns include at least time 
steps and the number of trial. The Boltzmann model 
analysis unit 5 receives the initial values and the evalu- 
ation conditions, and it repeats simulations of price fluc- 
tuation within the range of the specif '^ condition, based 
on the Boltzmann model using a Monte Cario method. 
A velocity/direction distribution setup unit 8 supplies the 
probability distributions of the price, the price change 
rate, and the price change direction for the financial 
product to the Boltzmann model analysis unit 5. The sys- 
tem also has a random number generator 9 for a Monte 
Cario method employed in the analysis by the Boltz- 
mann model, and an output unit 6 for displaying the anal- 
ysis result. 



A dealing system 100 applies the financial Boltz- 
mann model to option pricing, and reproduces the char- 
acteristics of Leptokurcity and Fat-tail by linear Boltz- 
mann equation in order to define risk-neutral and unique 
probability measures. Consequently, option prices can 
be evaluated in a risk-neutral and unique manner, taking 
into account Leptokurcity and Fat-tail of a price change 
distribution. 
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Description 

[FIELD OF THE INVENTION] 

[0001 ] The present Invention relates to a system for assessing a price distribution or a risk distribution for a financial 
product or its derivatives, which can rigorously evaluate a price distribution or a risk distribution, including a probability 
of occurrence of a big price change, based on a Bottzmann model. This system is also capable of analyzing price 
fluctuatk)n events for the financial product or its derivatives that could not be reproduced by the conventional technique. 
[0002] The present invention also relates to a dealing system used in the financial field. 

[0003] The present invention further relates to a computer-readable recording medium storing a price and risk as- 
sessing program for a financial product or its d^tvatives, and to a computer-readable recording medium storing a 
ctealing program. 

[TEGHNOLOGICAL BACKGROUND] 

[0004] A technique for analyzing past records of price change in a financial product or its derivatives and for sto- 
chastk:ally obtaining a price distributbn or a risk distribution is generally called a financial engineering technology 
[(K)05] In general, the Wiener process is used to model a change of stock price in the conventional financial engi- 
neering technology. The Wiener process is a type of the Markov process, whk:h is a stochastk: process on condition 
that a future state is independent of a past process. The Wiener process is often used to descn'be the Brownian motions 
of gas-molecules. 

[0006] With variables of t (tinne) and z depending on the Wiener process, the Wiener process is characterized in the 
following relatbnship between At and Az that is an infinitesimal change in z during the infinitesimal time A t 



here e is the random sample from the standard Gaussian distribution. 

[0007] Thus, the Wiener process evaluates fluc^ations with the variables based on the standard Gaussian distribu- 
tion. 

[0008] The conventional risk evaluatton nrtethod for a financial product or its derivatives generally establishes upon 
applying the Ito's process, whk:h is developed from the Wiener process. The Ito's process adds adrift temn to the Wiener 
process on the assumption that a change of stock prk^ follows the Wiener process, and further introduces a parameter 
function of time and other variables. 

[0009] The price change in stock priced expressed by the Ito's process is 



here S is the stock price, r is the non-risl^ interest rate, a is the volatility (i.e., the predk^ed change rate), and W is the 
normal distrbution with the expectation value of zero (0) and the standard deviation of one (1). 
[0010] The simplest example of the ito's process is the geometric Brownian motion model of stock prices. With the 
geometrk: Brownian motion, equation (2) becomes 



Az = e7At 



(1) 



dS = rSdl + So 75w 



(2) 




(3) 



here x is the natural logarithm of the stock price S. 

[0011] The probability density furK:tion Ptx; t) ol x based on equation (3) is 
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(4) 



[0012] Equation 4 is the Fokker-Plank equation, and is a typical diffusion problem. The solution of equation (4) is 



and the probability density function P(x; t) of x becomes the Gaussian distribution. 

[0013] Equation (5) is characterized by not only its simple form, but also effectiveness in evaluating price changes 
for financial products, because it is known that the price change of derivatives derived from the underiying assets has 
the same shape as that of the underlying assets (Ito's theorem). For this reason, various financial derivatives have 
25 been reproduced. 

[0014] However, the conventional technique for evaluating risks for the financial products or the derivatives is not 
capable of providing suffctentty reliable results, as is known in this field. 

[0015] This is because that conventional method evaluates risks of financial products based on the Gaussian distri- 
bution, and therefore, the probability of occun^ence of a big price change is underestimated. 
30 [0016] Although the likelihood of occun^nce of a big price change is low, such a big price change has a signifk^nt 
influence to investing risks, as compared with situations under the nomrial price changes. Accordingly, any risk evalu- 
ation methods or systenrts for financial products can not be reliable in the practical aspect unless the probability for the 
big price changes is accurately treated. 

[0017] Another problem is that the conventional risk evaluation technique requires some corrections to a heteroge- 
35 neous problem, in whk:h the probability density function changes depending on prices, or to a non-linear problem, in 
whbh the probability density functk>n used for the evaluation is a non-linear function. Along with the conventional 
approaches, such corrections have to be added empirically or based on know-how. In other wonis, the conventional 
technique requires dealer's experiences or uncertain judgements in the actual nrtaricet trading. 
[0018] Furtherrnore, the cofiventional risk evaluatk>n technique has; very limited capabiliti^ definition, 
40 and selection of the variables to produce price fluctuations observed in the markets. In other words, with the conven- 
tional technique, the probabOity density functbn can not be sufTidentty evaluated with variables for describing risks for 
financial products, if the actual price change distribution of an financial product is located out of the standard Gaussian- 
type distribution. This insuffk^iency can also be true in the cases where the price change mte is influenced by the past 
price change rate, and correlations exist b^ween the probability for price-up and the probability for price-down, or 
45 between the price change rate and the price change direction. The conventional technique is not capable of describing 
the probability density functk>n for the price change directk>n as well, and therefore, the probability distribution of the 
price change direction for the financial products are cfisregarded. 

[0019] Further problems in actual application of the conventional risk evaluation method for financial products relate 
to insufficient numerical techniques, such as random number sampling and variance reduction areas for making Monte 

50 cario cateulations. Consequently, undesirable variance inevitably remains in the conventional technique. 

[0020] Meanwhile, dealers or traders use dealing systems in banks or security companies for purposes of supporting 
the transactk>ns. The conventional dealing systems catoulate theoretical prices of financial products or its derivatives 
(hereinafter, referred to as "optrons"), and simulate risk evaluation or position change based on accepted theories, 
such as the Black-Sholes model or its expanded models. The Black-Sholes model assumes that the probability distri- 

55 bution of, for example, a stock price at a feature point of time is the Gaussian (or normal) distribution (Black, F. & M. 
Sholes, The Pricing of Optk>ns and Corporate Liabilities", Journal of Political Economy, 81 (May-June 1973), pp. 
637-59). 

[0021] However, the conventbnat dealing systems have many problems listed below in item (1 ) through (8). 
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[0022] (1) A so-called Fat-Tail problem Is a serious problem in the financial field (Alan Greenspan, "Financiat deriv- 
atives"; March 19. 1999; httpy/Www.federalreserve.gov/boarddocs/speeches/1 999/1 999031 9.hlm) . 
[0023] In order to calculate theoretical prices for financial products or options or to simulate risk evaluation and 
position change, the usage of the Gaussian distribution In the conventional financial engineering models has facilitated 

5 evolution of theoretical financial engineering and implementation to a computer system easy-applicable in business 
areas. In other words, non-normality of probability distributions have been often observed financial markets, in which 
big price changes actually occur, or in whbh transa<^ons are not so ac^ve. Once if this fact would be introduced from 
the early state, evolution theoretical financial engineering deployment and implementation to a computer system would 
have become much more diffbult in actual applk^ation. For this reason, dealers have to transact relying upon their own 

10 experiences or intuitions. To carry out such trartsactions, it is critk:alty important for the dealers to accurately grasp the 
behaviors for volatility of the market. 

[0024] (2) There ts-volatility defined as a historical volatility that is cateulated tyy the observed fluctuatk>n of prrce 
under asset transactions. A general method for obtaining a histork:al volatility is to calculate the standard deviation of 
returns for the asset, based on the observed fluctuation of the closing price. Other known methods for cabuiating a 
IS historical volatility include so-called "Extreme Value Theory" for estimating a volatility from daily high prices and low 
prices, and a modified Paricinson method for estimating a volatility while taking discontinuity of time in the actual trans- 
ac^'on Into account. 

[0025] However, these methocte have drawbacks in actual apprication for some reasons. For example, if transaction 
is not active enough, there is no continuity In the movement for the price of the underlying assets. In case that the 
20 dosing prices are used for determining the historical volatility, the corresponding exercise time wouki be altered from 
the actual one. 

[0026] Even in the case that the transaction is active, the conventional methods are not suitable to estimate the 
volatility under the fat-tailed regime because these methods assume nomnality in the risk probability distn*bution for the 
maricet behaviors. To this end, the volatility cabulated by the conventk>nal methods is used only as a rough guideline 

25 within the limited applk^ability. 

[0027] (3) Implied volatility (abbreviated as "IV*) is known in the option market, other than the historical volatility 
mentioned above. Implied volatility is volatility cak:ulated back from the option prices obsen^ed In the market along 
with the Black-Sholes equation. Implied volatility Is often used as a facto rfor calculating the theoretical price of an option. 
[0028] However, in a non-active market (for instance, the option transaction market for underiying assets in the 

30 ' current Japanese security market), as the number of observed transactions for options is small, the implied volatility 
for the conresponding options cannot be well-defined from the actual market data. For this reason, dealers have to 
repeatedly manipulate factors for regular adju^'ng the volatility parameters through their own experiences and expert- 
judgements in order to reflect the current market behavk>rs in their dealing activities. 

[0029] (4) Implied volatility of an individual stock option can provide some kind of important information to know what 
35 is the market reaction to the volatility of a particular stock for dealers, while the maricet reaction varies as transactions 
go on reflecting market circumstances. 

[0030] (5) In the option market, different sets of the implied volatility are obtained for multiple opttons originated from 
the same original underlying asset. This phenomenon is referred to as a "smile effect". I n this case, approximated smile 
curves are drawn in the two-dinnensionat phase space with the combinations between a vertical axis showing the 
4a magnitude of implied volatility and a horizontal axis showing the exercise prices for the options. Those smile curves 
are often used to cabulate theoretical prices of the option. 

[0031] A volatility matrix, whk:h is a data table having a time dimension along maturity of the option, in addition to 
the two-dimensional phase space nnentioned above, provkie infomriation to obtain the theoretical exercise price and 
to interpolate the volatility value for the regions unobserved in the market behaviors up to maturity, in conjunction with 
45 the above-mentioned Item (4). 

[0032] (6) In fact, in order to obtain reliable smile curves or a volatility matrix, suffk^ient nunnbers of option prices 
must be observed in the market. On the other hand, in some cases that under the transactions in moderately active 
market, the optk>n prices observed in the market are scattered in a wide range, it becomes difTicult to grasp a compre- 
hensive trend. 

50 [0033] (7) In general, the amount of information observed in the market trades off the rationality needed in the bask: 
assumption of the required model for estimating theoretk:a) price for the options relating to the underlying asset price. 
If numbers of the past records for the option prunes observed in the market are insuffk:ient, a stronger assumption is 
required in a model used to obtain dynamks of an expected probability distribution for the underlying asset price. As 
one of the advanced riKxiels having strong assumptions, the stochastk; volatility model (SVM) (Hull, John C. & Allan 

55 White, "The Pricing of Options on Assets with Stochastk: Volatilities". Journal of Finance, 42, June, 1 987, pp. 281 -300), 
and the GARCH model (T. Bollerslev, "Generalized Autoregressive Conditional Heteroskedasttctty". Journal of Econo- 
metrics, Vol.31 , 1986, pp. 307-327) are well-known. However, because these models assume nonmality in the proba- 
bility distributions, th^ can not deal with the Fat-Tail problem sufficiently. 
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[0034] There is a technique expanded from a lattice method and having no assunnption of normality (Rubinstein, 
Mark, "Implied Binomial Trees", Journal of Rnance, 49, July 1994, pp.771 -818). This technique is capable of forming 
a flexible probability distribution talcing the smite effect into account However, this technique requires sufficient numbers 
of option prices observed in the market in order to detenntne the distr&ution fomi. Therefore, this technique is not 

5 suitably used in a non-active option market 

[0035] There is also a famous model named a jump model, which independently generates a stochastic process 
entirely different from the nonnal distribution for a Fat-Tail problem (e.g., R.C. Merton, "Option Pricing When Undertying 
Stock Returns Are Discontinuous", Journal of Rnandal Economics, Vol.3, March 1976, pp.1 25-1 44). However, the 
Jump model has an assumption of discontinuous price fluctuation, and therefore, the stochastic volatility model (SVM) 

10 naturally becomes a nonlinear prot>lem. For this reason, the risk-neutral measure can not be achieved invariably, which 
prevent the option price from being d^ned uniquely. 

[0036] (8) In conduskm, no conventional technique can provide minute and accurate infonnation in real time for 
solving the Fat-Tail problem and being applicable to a non-active market, although it has been desired for dealers and 
traders to receive significant smile curves or a volatiltty nriatrix on their displays in real time in response to the actual 
IS market that changes every nrK>nDent. The conventional technique is incapable of automatically acquiring necessary 
data required for computation in response to requests from the dealers in an interactive manner, and of automatk:ally 
selecting the optimum model to analyze the nnarket deeply and ftexbly. 

[0037] Therefore, an object of the present inventbn is to introduce a probability density function with a higher accuracy 
in comparison with the nonnal distribution, and to develop a system capable of correctly evaluating the price distribution 
20 and the risk distribution for a financial product or its derivatives. 

[0038] Another object of the present invention is to provide a price and risk evaluation system for a financial product 
or its derivatives, which system is capable of theoretkrally sohnng the above-mentioned heterogeneous or nonlinear 
problems. 

[0039] It is still another object of the present invention to provide a price and risk evaluation system for a financial 
25 product or its derivatives, which system introduces a new function of probability density for estimating a price distribution 
and a risk distribution. This probability dertsity function model can adequately define and describe variables that can 
not be dealt sufficiently by the conventional technique, and can establish a reliable method. 

[0040] It is still another object of the present invention to provide a price and risk evaluation system for a finandal 
product or its derivatives, which system Introduces a new function of probability density for estimating a price distribution 
30 and a risk distribution of a finandal product This function is capable of establishing a sampling method for improving 
the eftidency of computation, and alfov^ risk prk^es to be connputed at a high effidency. 

[0041] It is yet another object of the present invention to provide a price and risk evaluation system for a finandal 

product or its derivatives, which system ts applicable to a parallel computing with high efficiency. 

[0042] It is yet another object of the present invention to provide a computer-readable recording medium storing a 

35 dealing program, which includes a Bottzmann model computation developed by the nuclear reactor theories and applied 
to the finandal field, in place of the general theories used in the conventional techniques. This program is capable of 
dealing with big price changes in the undertying assets (a fat-tail problem mentioned above), and is applicable to an 
option market in which transactions are not so active. This program allows a computer system to display significant 
theoretk:al prices arKi risk parameters on display terminals of dealers and traders by means of the Interactive screen 

40 interfaces. 

[SUMMARY OF THE INVEISTTION] 

[0043] The first object of the present invention is to provkie a price and risk evaluation system for evaluating a price 
^ distribution or a risk distribution for a finandal product or its derivatives with applications of a new calculation method, 
the details of whk^h will be described bek>w. A risk cak^ulation method is also provided, which is capable of treating the 
fat-tail problems in the finandal engineering field assodated with underiying assets and their derivatives, and of sub- 
stantially eliminating defects or drawbacks having existed in the prior art. 

[0044] This system has an initial value setter and an evaluation condition setter. The initial value setter receives at 
so least one of the initial values of a price, a price change rate, and a price change direction for a financial product or its 
derivatives that be evaluated. The evaluation condition setter allows a user to input evaluation conditions including at 
least one set of tinr^ steps and the numt)^ of trials for cak:utations. As a significant feature, the system has a Bottzmann 
model analyzer, which receives at least one of the initial values and the evaluation conditions, and repeats simulattons 
of price fluctuation based on a Boltznnann model using a Monte Cario method within the ranges of the given cak:ulation 
55 conditions. The Boltzmann nnodel analyzer can obtain a price distribution or a risk distribution for the financial product 
or its derivatives in an accurate manner. The risk evaluation system also has a velocity/direction distribution setter that 
supplies probability distributions of the price, the price change rate, and the prrce-change direction for the finandal 
product or its derivatives to the Bottzmann model analyzer. This system has a random nunnt>er generator in the Boltz- 
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mann model analyzer, and an output unrt that outputs series of analysis results from the Boltzmann model analyzer. 
[0045] The initial value setter acquires the initial valu^ of the price, the pnce change rate, and the price-change 
direction for the financial product or its derivatives from a mari(et database that stores information about detailed trans- 
action histories, such as exercises and ask-bit data. The initial value setter then supplies the acquired initial values to 
5 the Boltzmann model analyzer. The velocity/direction distribution setter receives the past records of a selected financial 
product or its derivatives from the market database, and generates a probability density function with variables of the 
price, the price change rate, the price change direction, and time. The velocity/direction distribution setter then supplies 
the probability density function to the Boltzmann model analyzer. 

[0046] The price and risk evaluation system further has a total cross-section/ stochastic process setter, which supplies 
10 information for setting a sampling time width of the simulation of price fluctuation to the Boltzmann model analyzer. In 
this case, the total cross-section/stochastk:, process setter acquires a price fluctuation frequency and a price change 
rate of the financial product or its derivatives from the market database storing infomnatbn about financial products or 
derivative products. The total cross-section stochastic process setter then inputs a ratio of the price fluctuation fre- 
quency to the price change rate into the total cross-section term of the Boltzmann's equation. 
15 [0047] The velocity/direction distribution setter acquires the past records of a selected finandal product or its deriv- 
atives from the market database storing infonnation. The vek)city^direction distnl)ution setter than infers a distribution 
of the price change rate forthefinanctal product or its derivative using a Sigmoki function and its approximation fonns, 
and supplies the infen'ed distribution of the price change rate to the Boltzmann model analyzer. 
[0048] Alternatively, the velocity/direction distribution setter acquires the past records of a selected financial product 
20 or its derivatives from the market database, and determines a set of the S^otd function parameters to estimate the 
probability distribution of the price change rate using the price change rate data stored in the past records. The vekx^/ 
direction distribution setter then supplies the distribution to the Boltzmann model analyzer. 

[0049] In still another alternative, the velocity/direction distribution setter acquires the past records of a financial 
product or its derivatives from the market database, and estimates a probability distribution of the price change direction 
25 for the financial product or its derivatives using the past records. The velocity/direction distribution setter then supplies 
the probability distribution of the price change direction to the Boltzmann model analyzer. 

[0050] This vek>city/direction distribution setter infers the probability distn'bution of the price change direction, taking 
into account a correlation behween the probability for price-up and the probability for price-down. 
[0051] In still another altemative, the velocity/direction distribution setter acquires the past records of a financial 
30 product or its derivatives from the maricet database, and generates a probability distribution of the price change direc- 
tion, taking Into account a correlation between the distribution of the price change rate and the distribution of the price 
change direction for the financial product or its derivatives. The velocity/direction distribution setter then supplies the 
probability distributions to the Boltzmann model analyzer. 

[0052] In still another altemative, the velocfty/direction distribution setter generates homogeneous probability distri- 
35 butions ttiat are independent of the prices, or heterogeneous probability distributions that depend on the prices, with 
regard to the probability of a price change rate and a price change direction distributions. The velocity/direction distri- 
bution setter then supplies these homogeneous or heterogeneous probability distributions to the Boltzmann model 
analyzer. 

[0O53] In this system, the Boltzmann model analyzer obtains the price distn'bution or the risk dstribution for the 
^ financial product or its derivatives using either a linear Boltzmann model or a non-linear Boltzmann model. In the linear 
Boltzmann model, the cross-section used in the Boltzmann's equation is independent of a probability density or flux 
for the financial product or its derivatives. In the non-linear Boltzmann nrmdel, the cross-section for the Boltzmann's 
equation is dependent on the prot>abilrty density or the flux for the financial product or its derivatives. 
[0054] Altemativety, the Boltzmann model analyzer obtains the price distribution orthe risk distribution for the financial 
^ product or its derivatives using a product of a probability density function and a price change rate per unit time for the 
financial product or its derivatives, as flux of the Boltzmann's equation. 

[0055] The Boltzmann model analyzer evaluates a probability density at an arbitrary time based on the track-length 
calculated using flux for the finartctal product or its derivatives in order to reduce a variance. 
[0056] The Boltzmann model analyzer evaluates a price probability in an infinitesimal price-band or a risk probability 
so in an infinitesimal time inters using all of or a part of the price fluctuation data f orthe financial product or the derivatives. 
Thus, the Boltzmann model analyzer reducesa variance of the price orthe risk by applying the point detector technique 
often empbyed in a neutron transport Monte Cario simulation. 

[0057] The Boltzmann model analyzer cak:ulates an adjoint probability density or an adjoint flux deduced from an 
adjoint Boltzmann equation for a price fluctuation of the financial product or the derivatives, and reduces variance by 
55 weight-sampling technk^ue using values in proportion to the adjoint probability density or the adjoint flux. 

[0058] In the price and risk evaluatton system for a financial product and its derivatives, the velocity/direction distri- 
bution setter generates a vebdty distribution or a direction distiribution for a financial product or its derivatives, taking 
into account the inter-correlation among other multiple financial products or their derivatives. The generated probability 
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distribution functions are supplied to the Bottzmann anatyzer. 

[0059] The Boltzmann analyzer evaluates a price distribution or a risk distribution for a financial product, and then, 
applies the Ito's theorem to calcu late the equivalent price distribution or the equivalent risk distribution for the derivative 
of that financial product. 

5 [OOGO] Preferably, the Boltzmann analyzer consists of multi-methods for carrying out simulation of price fluctuations 
and constructs the edited probability density function through gathering the simulated price fluctuations. 
[0061 ] According to the present invention, a price or a risk evaluation system of the first embodiment described above 
can correctly evaluate the probability of occunrence of big price changes for a financial product, compared with the 
conventional technique applying the nomDai distribution as for the probability density. Consequently, the price distribu- 

10 tion or the risk distributk>n for the financial product or its derivatives can be evaluated more accurately 

[0062] The present system can treat probability densities in heterogeneous problems in which the probability density 
functk)n used for evaluation varies depending on the price, or in nonlinear problems in which the probability density 
functk)n itself is nonlinear, without heavily relying on experiences or know-how. 

[0063] The present system also can evaluate probability distributions in more flexible man ners than the conventional 
IS methods, especially in the following situations: 

non-Gaussian probability distributions, which were not property taken into account in the conventional methodol- 
ogy, are treated; 

price probability distributions are influenced by the past price itself; 
20 price probability distributions are characterized by the inter-correlation between the probabilities of price-up and 

price-down; and 

probability distributions are conflated with the price change rates and the price change direction. 

[0064] The present system requires no time grid for simulatk>n of prk:e fluc^ation, whereas the conventional tech- 
25 niques need time grids for the calculations. The system can evaluate a probability distribution at an arbitrary point of 
time within the observed area by introducing the flux concept, whereas the conventional technique can treat a probability 
distribution only at a selected time. 

[0065] The present system ateo introduces an idea similar to the point detector used in the neutron transport calcu- 
lation by Monte Cario simulations. The Idea of point detector allows the system to automatically detect a route of 
30 causing an event in a target region in an infinitesimal observation area, in which no price change occurs or no flux 
passes, based on all of or a part of the price-changing events. Accordingly, the system can evaluate an even in an 
arbitrary infinitesimal observation region, while reducing the variance. 

[0066] The present system also reduces the variance generated in the Monte Cario cabulation for a probability 
density by introducing the idea of adjoint probability density or adjoint flux into the financial technology and selecting 

35 weights in proportion to the magnitudes of adjoint flux in the phase space. 

[0067] In the situatbns where inter-correlations are found among the probability distributrons of multiple financial 
products, the present system can evaluate a price distribution or a rek distribution taking such con^lattons into account 
The system includes applk:ation of the lto*s theorem to evaluate the probability distribution of the price change rate for 
derivatives based on the probability distribution for the underlying asset 

40 [0068] The present system is realized with parallel computatk>n systems. Consequently, a price distribution or a risk 
distribution for a financial product or its derivatives can be obtained in a highly efficient manner based on the parallel 
computation. 

[0069] The present system introduces the Boltzmann model, instead of diffusion models having an assumption of 
the standard nonrnal distribution used in the conventional systenrts. Accordingly, the present system can be replaceable 
^ or substitutable for the existing financtal-retating systems for evaluating risks or analyzing portfolio. This means that 
the hardware resources and the various types of infonmation required in the existing system and installed by the con- 
ventional methodology can be utilized as they are. As a result, an efficient system for evaluating a price and a risk 
probability distribution for a financial product or its derivatives can be realized. 

[0070] As the second embodiment of the present invention, a computer-readable recording medium storing a program 
50 of prbe and risk evaluation is provided. By installing this program into a computer and causing the computer to execute 
the following operatk>ns, a prk^e/risk evaluation system can be built up. Namely, an initial value setter of the computer 
system inputs at least one of the initial values of a price, a price change rate, and a price change direction for a financial 
product or its derivatives. An evaluation conditk>n setter of the connputer inputs evaluation conditions including at least 
time steps and the number of trials. A Boltzmann model analyzer of the computer repeats simulations of price fluctuation 
55 using a Monte Cario method based on the Bottzmann model within the ranges of the evaluation conditions in order to 
obtain a price distributton or a risk distribution for the financial product or its derivatives. A velocity/direction distribution 
setter supplies the probability distributions of the price, the price change rate, and the price change direction for the 
financial product or its derivatives to the Boltzmann model analyzer. A random number generator yields a series of 
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random numbers used for a Monte Carlo analysis in the Bottzmann model analyzer, and an output unit provides various 
types of outputs from the analysis result obtained by the Bottzmann analyzer. 

[0071] The third embodiment of the present invention covers a computer dealing system applying the Bottzmann 
method and its related calculation tools. The dealing system comprises an implied volatility computation engine, a 

5 Bottzmann model computation engine, a conversion fitter, and a dealing temnrnal. The implied volatility computation 
engine provides an implied volatility based on market data. The Bottzmann model computation engine evaluates an 
option price for a selected option product based on the Bottzmann model using the market data. The conversion filter 
of the dealing systm converts the option price obtained by the Bottzmann model computation engine into a volatility 
of the Black-Sholes equation. The dealing terminal provides displaying the volatility of the Bl^k-Sholes equation in 

10 comparison with the implied volatitrty cak:utated from the market data, or displaying the option price calculated by the 
Bottzmann model computation engine in comparison with an option price in mari(et. 

[0072] The present dealing system define a unique and risk-neutral probability measure by applying the Bottzmann 
model used in the financtal engineering field to option price evaluation, because the system can treat Leptokurcity and 
Fat-tail problems for a price or a risk probability distrft)ution appropriately in a linear equation fomn. Consequently, the 
IS system can evaluate option prices with the risk-neutral and unique manner taking into account the Leptokurcity and 
Fat-tail of the price-change distribution. Applying the Bottzmann model to the option price evaluation of a selected 
option product allows the system to grasp the comprehensive tendency of a volatility matrix from the past transaction 
records varying in a wide range. 

[0073] The present dealing system covers dealing of a stock price index option or an IndMdual stock option as the 
20 option product Consequently, the comprehensive tendency of the volatility matrix of the individual stock option can be 
obtained. 

[0074] As described above, the present dealing system provides whole trends for volatility matrices for the individual 
stock option whose dealings are less active, by checking the consistency of the daily eaming rates with the conre- 
sponding underlying assets. This can be achieved because the Bottzmann model is capable of pricing the optk)ns 
25 through detennining a set of Bottzmann parameters as to reproducing the daily eaming probability reflecting the market 
data for the underlying assets. 

[0075] The Bottzman n model computatk>n engine of the present dealing system has a cateulation unit that cabulates 
an optk>n price consistent with historical infonmatbn. Accordingly, a well-adjusted option price can be provided to the 
user via the dealing terminal. 

30 [0076] In the present system, the Bottzmann model computation engine also has a converter that converts option 
prices, whk^h were sets of exercise prices for the discrete months of the delivery, into sets of equivalent volatility from 
the Black-Sholes equation. The equivalent option prices and the risk parameters are obtained through interpolation of 
the Black-Sholes equation, and are displayed on the dealing terminal of the user 

[0077] The Boltzmann model computation engine In the present system has a table generator that generates a table 
35 of a prot>abitity density function evaluated by the Boltzmann model, and calculates an option price from the sum of 
inner product of vectors (i.e., Riemann sum). This arrangement increases the operation speed, and realizes a highly 
responsive ctealing system. 

[0078] As the fourth embodiment of the present invention, a dealing system comprising a dealing terminal, a theo- 
retcat option price and paranrteter computation engine, an interpolatk>n unit, and an interface that receives market 

"fo data. The dealing tenminal functk>ns as a graphics user interface. The theoretk^al option price/parameter conoputation 
engine is configured to switch between a rough computatk>n and a detailed conrtputation of a theoretical option price 
and parameters. The rough computation provides the theoretk:al option price and the parameters for each exercise 
price and for each delivery month set in the normal market-state. The detailed computation provides more detailed 
information including the theoretk^al option prices and the parameters for those exercise prices and the specified de- 

^ livery months by users that are not set in the market, but designated by a user. The price interpolation unit covers 
pricing with the Bottzmann model up to artntrary maturities specified by users. This dealing system provides the rough 
. computation results to the dealing temninal as a market activity at a high speed in the ordinary state. The dealing system 
displays the detailed computatk>n results on the dealing terminal for a specific price band ^edfied by the user 
[0079] In the present dealing system, quk* martcet display is made for the rough computation results in the normal 

50 state, rather than displaying the detailed computation results. On the other hand, the detailed computation results are 
displayed to timely inform the user of any changes in the maricet activity, by calculating the detailed theoretical price 
and parameters based on the Bottzmann model within the associated price range in response to the user's instruction. 
[0080] As an alternative, the dealing system conrtprises a dealing tenminal. a rough computation engine, a multi^erm 
Boltzmann engine, an interpolation unit, and an interface. The dealing terminal functions as a graphics user interface. 

55 The rough computatk>n engine computes a theoretcal option price and parameters for an exercise value of each 
delivery ntonth set in the market The nrujlti-temn Boltzmann engine computes theoretfcal optton prices and parameters 
at arbitrary terms based on the Bottznnann nruxtel. The dealing system nonnalty causes the dealing terminal to display 
the market activity based on the rough computation result, and causes the display temninal to show the multi-term 
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volatility in response to the user's instruction. The rough computation results are nnade faster than the detailed com- 
putation results. Because the present system can evaluate the term structure of volatility that does not come up in the 
market data until the end of time concerned, the developing efficiency of a structured bond or an exotic option can be 
improved. 

s [OOai] As the fifth emt>odiment, a dealing system, which comprises a dealing temiinal having a graphics user inter- 
face, a rough computation engine, a detailed computation engine, an interpolation unit, a position setter, an automatic 
transaction order unit, and an interface for receiving market data, is provided. The rough computation engine computes 
a theoretical price and an index for each exercise price and for each delivery month set in the market. The detailed 
computation engine computes detailed irrfomiation including theoretical pnces and parameters for exercise pnces and 

'0 ctelivery months that are not set in the market This dealing system outputs an automatic order signal when a stock 
index option price or an individual stock option price reaches a predetemnined autonnatic ordering price band. TTiis 
system allows the user to visually confirm the appropriate standard with an advanced model, to set a position, and to 
timely order in an automatk: manner. 

[0082] The present dealing system with the graphical user interface (i.e., the dealing terminal) facilitates a fading 
'5 processor as an altemative. Inthtscase, the dealing system causes the dealing terminal to display an animated behavior 
with a fading style for a term structure of a volatility that has been converted from an option price at the money (ATM) 
obtained from the Boltzmann model. This an^angement prevent mtspridng due to overreactions to the maricet 
behaviors. 

[0083] The present dealing system installs a risk limit setter. In this case, the dealing system alerts the user when 
20 the price in the market enters the warning area specified by the user. In other words, the system allows the user to set 
a risk limit, and ermbles the dealers cortcemed to conduct risk management appropriately. 

[0084] The present dealing system installs an alternative-position selector, in addition to the risk limit setter unit. In 
this case, the dealing system alerts the user, and simultaneously causes the dealing terminal to display an alternative 
position. Since the alternative position is automatically selected, the system can prevent the user from suffering a loss 

25 due to overreactions to the fluctuations in the market price. 

[0085] As the sixth embodiment of the present invention, a corinputer-readable recording medium storing a dealing 
program is provided. By installing this program in a computer and causing the computer to execute the following op- 
orations, a dealing system can be built up. Such a dealing system computes an innplied volatility and an option price 
of a selected asset based on the Boltzmann model using the market data. The system then converts the option price 

30 obtained from the Boltzmann model into an equivalent volatility from the Btack-Shples equation. Finally, the system 
displays the equivalent volatility from the Btack-Sholes equation in comparison with the implied volatility calculated 
from the nnarket data, or to display the option price based on the Boltzmann model in comparison with an option price 
in market 

35 [BRIEF DESCRIPTION OF THE DRAWINGS] 
[0086] 

Fig. 1 represents a block diagram showing the structure and the operation flow of a price and risk evaluation system 
"10 for financial products according to the present inventk)n; 

Rg. 2 illustrates an operation flow of the Boltzmann model analysis unit according to the present invention; 
Fig. 3 illustrates another operation flow of the Boltzmann rruxlel analysis unit which uses a price chemge frequency; 
according to the present invention; 

Fig. 4 scherriaticalty illustrates the simulation results of the operation flow of Fig. 2 in the observation areas; 
^ Rg. 5 schematk:alty illustrates the simulation results of the operation flow of Fig. 3 in the observation area; 

Rg. 6 is a graph showing a profc>ability distribution simulating a diffusion model using the Boltzmann model accord- 
ing to the present inventk>n; 

Rg. 7 is a graph showing price changes simulating the diffusion model using the Boltzmann model according to 
the present invention; 

50 Fig. 8 is a graph showing required spectra with respect to the price change rate v of stock prices; 

Rg. 9 ^ a graph showing the dependency of the spectra on the inckient velocity v*; 

Fig. 10 is a graph showing the dependency of the price-up component (positive changes) of the spectra on the 
incttlent vek)Cfty v*; 

Rg. 1 1 is a graph showing the dependency of the price-down component (negative changes) of the spectra on the 
55 incident velocity v*; 

Rg. 12 is a graph showing an application of evaporation spectra; 

Rg. 13 is a graph showing the empirical equatton of the velocity term of the differential cross-section; 

Fig. 14 illustrates the averages of the continued price-up probability and the continued price-down probability of 
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every 5 days; 

Rg. 15 is a graph of the simutation result of price fliK^ation using a Boitzmann model; 

Fig. 16 shows the simulation resuft of price fluctuation using the Boitzmann model with a detailed view; 

Fig. 1 7 illustrates the distribution of the stock price after two hundred days using the Boitzmann model; 

Fig. 18 is a graph of stock price distributions of every twenty days using the Boitzmann model; 

Fig. 19 illustrates the configuration of parallel processing of the price and risk evaluation system for a financial 

product or its derivatives according to the present invention; 

Fig. 20 is a graph showing the price fluctuation CI of the underlying assets expected by the geometric Brownian 

model, in comparison with the ck>sing-price fluctuation (daily earning rate) C2 of a typical stock price; 

Rg. 21 is a graph of fluctuation of the daily earning rate C3 of the stock price average of the Nikkei 225 stock 

average; 

Rg. 22 is a graph showing the implied volatility of the put option of the price index of stocks of the Nikkei 225 stock 
average, together with a smile curve; 

Fig. 23 is a graph showing the probability density of an actual daily eaming rate, together with the normal distribution 
presunrted by Black-Sholes equation; 

Rg. 24 s a graph showing the price change probability estimated by the Boitzmann model, together with the 
logarithmk: nomriat distribution used in Black-Sholes equation; 

Rg. ^ is a graph of temperature T as a function of daily eaming rate v* of the previous day; 

Rg. 26 is a graph of simulated probability density as a function of daily eaming rate for the price valuation using 

the Boftzmann model; 

Fig. 27 is a graph showing the inr^lied volatility of the Boitzmann model, in comparison with a jump model; 
Fig. 28 is a block diagram for the dealing system according to the present invention; 

Fig. 29 is a block diagram showing the operation flow and the structure of the Boitzmann model computation engine 
used in the dealing system shown in Rg. 28; 

Fig. 30 is a flowchart of theoretfcal computation carried out by the dealing system shown in Fig. 28; 

Fig. 31 is a graph of the temperatures T of various listed stocks of Tokyo stock exchange obtained by the dealing 

system of the invention; 

Fig. 32 shows examples of evaluatk>n for call option pricing of individual stock options, which is expressed by ratio 

of call option price to underlying assets as a function of ratio of exercise price to underlying assets; 

Fig. 33 shows examples of evEiluation of put option pricing of individual stock options, whbh is expressed by ratio 

of call option price to underlying assets as a function of ratio of exerdse price to underlying assets; 

Fig. 34 is a graph indicating the relationsh^ between the ratio of exercise price of a call option to underlying assets 

and the implied volatility, whbh is obtained by the dealing system of the invention; 

Fig. 35 is a graph indicating the relationship between the ratio of exercise price of a put option to underlying assets 
and the implied volatility, whch is obtained by the dealing system of the Invention; 

Fig. 36 shows the operation flow of theoretical computation flow carried out by the dealing system according to 
the invention; 

Fig. 37 illustrates a sub-screen displayed on a terminal of the dealing terminal of system, whk:h displays the detailed 
track of stock index in the continuous session; 

Fig. 38 illustrates a sub-screen displayed on a terminal of the dealing system, which displays a table of the implied 
volatility, themarket prices for each delivery month, and each exercise price of stock index option, together with 
the underlying assets; 

Fig. 39 illustrates sub-screens of the dealing system, whk:h displays, in graphs, information contained in the table, 

such as the implied volatility, the martcet prices of stock index options for each exercise price and each delivery 

month, using the stock indrces displayed on the terminal as an undertying asset price; 

Fig. 40 illustrates an operation flow of the detailed price evaluation canied out by the dealing system; 

Fig. 41 illustrates switching of graphs displayed on a terminal of the dealing system during the process of detailed 

price evaluation shown in Rg. 40; 

Rg. 42 illustrates switching of tables displayed on a tenrtinat of the dealing system during the process of detailed 
price evaluation; 

Fig. 43 is a graph representing the theoretical computation carried out by the Boitzmann model computation engine 
of the dealing system; 

Rg. 44 Illustrates an operation fk>w of the evaluation process of an arbitrary multi-temn volatility carried out by the 
dealing system; 

Fig. 45 illustrates a period setting screen displayed on a terminal of the dealing system in the evaluation process 
of the arbitrary multi-teon volatility shown in Rg. 44; 

Fig. 46 illustrates a table with market prk^es and innplied volatility values for the arbitrary multi-tenn delivery in the 
dealing system; 
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Fig. 47 is a graph showing an example of estimation result for the arbitrary multi-term volatility, which is displayed 
on the terminal of the dealing system as a result of the evaluation process; 

Rg. 48 illustrates an operation flow in a modified process forthe detailed price evaluation with a fading operation 
canied out by the dealing system, which includes a fading operation; 
5 Rg. 49 illustrates an operation flow showing the detailed steps for the fading operation; 

Rg. 50 is a table of the implied volatility and the option prices, in which a virtual implied volatility for a real-time 
ATM is inserted in the exercise price band set in the marked, used in the detailed price evaluation process canied 
out by the dealing system; 

Rg. 51 is a graph of implied volatility for arbitrary multi-terms, which includes the implied volatility for a virtual real- 
10 time ATM inserted in the exercise price band set in the market, displayed in the detailed price evaluation process 

carried out by the dealing system; 

Rg. 52 B a process diagram for explaining the fading operation; 

Rg. 53 illustrates an operatton flow in the arbitrary mutti-tenn volatility evaluation process using the fading function; 
Rg. 54 illustrates an operation flow for an automatic ordering process canied by the dealing system, in whk:h a 
'5 deal sets a desired position and orders timely; 

Fig. 55 is a graph showing the relationship between the output from the computation engine and the innplied vol- 
atility, which is used for setting a position; 

Rg. 56 illustrates a Input screen displayed on a tenntnal of the dealing system with the operation for position setting; 
Fig. 57 illustrates an operation flow in the flrst part of the process for displaying the simulating animation for the 
20 temn structure analysis of the ATM implied volatility in a fading manner in the dealing system; 

Rg. 58 illustrates an operation flow in the second part of the process for displaying the simulating animation for 
the term structure analysis of the ATM Implied volatility in a fading manner; 

Rg. 59 illustrates an operation flow in the third part of the process for displaying the simulating animation for the 
term structure aneitysis of the ATM implied volatility in a fading manner; 
25 Rg. 60 illustrates an operation flow for the detailed steps of the interpolation A shown in step S30 of Rg. 58; 

Fig. 61 illustrates an operation flow forthe detailed steps of the interpolation B show in step 830* of Rg. 58; 
Rg. 62 Is a graph showing the approach for setting the automatic ordering process; 

Fig. 63 illustrates an operation flow for setting an automatic waming func^on with market data in the dealing system; 
Fig. 64 is a graph showing the approach for setting the automatic warning process; 
30 Rg. 65 Illustrates an operation flow for automatically computing and selecting an altemative position, which is 

carried out together with the automate warning process in the dealing system; and 

Fig. 66 is a graph showlrtg the approach for setting the automatic waming process, together with selection of an 
altemative position In the dealing system. 

35 [Best Mode for Carrying out the Invention] 

[0087] The preferred embodiment of a price and risk evaluation system for a financial product or its derivatives will 
now be explained with reference to the drawings. Rg. 1 illustrates both the structure and the operation flow of the price 
and risk evaluation system for financial products. 

^ [0088] The price and risk evaluation system 1 1ncludes a portfolio setup unit 2, an Initial value setup unit 3 for inputting 
a price, a price change rate, and a price change direction, and an evaluation condition setup unit 4, a Boltzmann model 
analysis unit 5, an output unit 6, an total cross-section/stochastk: process input unit 7, a velocity distribution/direction 
distoibution setup unit 8, a random number generator 9, a V^R-evaluation unit 10, and a maricet database 11 . 
[0089] The Boltzmann model analysis unit 5 includes an initializatk>n unit 1 2, an initial value setting unit 1 3, a sampling 

^ unit 14, a price fluctuation simulation unit 15 for simulating price fluctu^ton baised on the Boltzmann model, a probability 
density calculating unit 16, a one-tria)-completk>n detector 17, an all-trial-completion detector 18, and a probability 
density editor 19. 

[0090] The rectangle defined by the long dashed line indcates the evaluation system 1 of the embodiment. The 
maricet database 11 and the VaR-evaluation unit 10 are positioned across the long dashed, which means that these 
so elements can be connected as external units to the evaluation system 1 via data communication. 

[0091] The long (tohed line does not intend to define a physk:al boarder of computer. For example, the elements 
included in the evaluation system 1 can be appropriately divided if the system is designed to cany out dispersed 
operations, such as a client-server system. 

[0092] The portfolio input unit 2 receives portfolio, and outputs a financial product or its derivatives for evaluation. 
S5 [0093] In an ordinary money management, the assets to be invested are allocated to multiple financial products or 
their derivatives In order to reduce a risk, and to carry out the most advantageous money management as a whole. A 
set of such multiple financial products or the derivatives, or the combination of these, is nanned portfolio. The portfolio 
input unit 2 extracts the financial product or its derivatives that is to t>e evaluated among fi^om the portfolio, and outputs 
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the extracted product. 

[0094] Preferably, the portfolio input unit 2 has a portfolio table or database Inside it, and allows a user to input the 
ID code of a desired portfolio. The portfolio setup unit 3 then exhibits the configuration of the portfolio, and allows the 
user to select a financial product or its derivatives concerned. 
5 [0095] The portfolio input unit 2 is not an essential element of the present invention, and therefore, rt may be omitted 
if the data required for evaluating the selected financial product or its derivative is known. 

[0096] The initial value setup unit 3 supplies at least one of the initial values of the price, the price change rate, and 
the price change direction for the financial product or its derivatives to be evaluated, to the Bottzmann model analysis 
units. 

10 [0097] The initial values of the price, the price change rate, and the price change direction for the financial product 
or its derivatives are obtained from the past record. Preferably, the initial value setup unit 3 receives the financial 
product or the derivatives from the portfolio setup unit 2, and retrieves information as to the financial product or its 
derivatives from the market database 11 . The initial value setup unit 3 then acquires the initial values of the price, the 
price change rate, and the price change direction for this financial product or its derivatives from the past records 

'5 contained in the retrieved information. The acquired initial values are supplied to the Boltzmann model analysis unit 
5. The initial-value input unit 3 is an essential element of the price/risk evaluation system 1 . 
[0098] The evaluatbn condition setup unit 4 supplies the evaluation conditions to the Boltzmann model analysis unit 
5. The evaluation conditk)ns include, for example, the number of trials, the time zone, and the price band for evaluation, 
whfeh are required for analysis by the Boltzmann model analysis unit 5. The evaluation condition setup unit 4 is an 

20 essential element of the evaluation system 1 , and supplies any evaluatbn conditions required by the Boltzmann model 
analysis unit 5. 

[0099] The Boltzmann model analysis unit 5 is the center of the price and risk evaluation system 1 . The Boltzmann 
model analysis unit 5 receives the initial values and the evaluation conditions from the initial value setup unit 3 and the 
evaluation condition setup unit 4, respectively, and repeats price fluctuation simulations for the selected finandal prod- 
25 ucX or its derivatives. The price fluctuation simulation is carried out based on the Boltznrtann model within the range of 
evaluation condition, using the Monte Gario method. The Monte Cario method is a numerical analysis method for 
obtaining a rigorous solution of a Boltzmann equation. 

[0100] The initialization unit 12 of the Bottzmann model analysis unit 5 initializes the values of the price, the price 
change rate, and the price change directk>n for the financial product and its derivatives to starting evaluation. 
30 [0101] The initial value setting unit 13 of the Boltznnann model analysis unit 5 sets up the initial values of the price, 
the price change rate, and the price change directk)n of the financial product based on the inputs from the initial value 
setup unit 3. 

[0102] The sampling unit 14 of the Boltzmann model analysis unit 5 detemiines a sampling widtti of the price fluc- 
tuation simulation. As the feature of the present invention, the sampling unit 14 can set a probability of price change 
35 per unit time based on the input from the total cross-section/stochastfc process setup unit 7. This arrangement allows 
setting of a time grid for simulation to be omitted. Setting of a time grid was required for simulation in the conventional 
system, but is a diffcult technique. This will be described in more detail below. 

[0103] The price fluctuation simulation unit 1 5 of the Boltzmann model analysis unit 5 simulates the price distribution 
from the last price based on the probability distributions of the price change rate and the price change direction, using 

40 the Monte Carlo m^od. 

[0104] The price fluctuation simulation unit 1 5 receives the price change rate and the price change direction for the 
flnarx^ial product from the price change rate distribution and the price change direction distribution setup unit 8, in order 
to simulate a price distribution based on a Boltzmann model. The price change rate and the price change direction 
conr^pond to the velocfty/direction distrftnitions in the Bottzmann equation. 

43 [0105] The prce fluctuation simulation unit 15 also receives a series of random numbers generated by the random 
number generator 9, and computes the solution of the Boltzmann equation by the Monte Cario method. 
[0106] The probabinty density cak:ulating unit 16 of the Boltzrr^n model analysis unit 5 integrates the price distri- 
bution simulated by the price fluctuation simulation unit 15 to obtain the probability density. 

[0107] The one-trial completion detector 1 7 of the Bottzmann model analysis unit 5 determines whether a trial has 
50 been completed. In this context, "one triaP means a series of simulation for price fluctuations in the time period from 
the beginning to the end of the simulation process. The one^trial completion detector 1 7 can detect the completion of 
one trial by comparing the curr^t point of time in the process simulation. The condition for completion of one trial is 
supplied from the evaluation condition setup unit 4. 

[0108] If one trial has not been connpleted yet for a series of simulation, the process returns to the sampling unit 1 4 
55 to caknjlate the price and the probability density for the subsequent time regime using the previous price, and distri- 
butions of the price change rate and the price change direction. 

[0109] The all-trial completion detector 18 of the Boltzrr^n model analysis unit 5 detemnines whether or not the 
current nunrtber of trials has reached the selected numt>er of trials given by the evaluation condition setup unit 4. The 
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maximum triaJ number is supplied to the all-trial completion detector 1 8 from the evaluation condition setup unit 4. 
[0110] The probability density editor 19 of the Bottzmann model analysis unit 5 collects the probability densities 
obtained from the entire trials and edits the probability density of the price fluctuation for the financial product or its 
derivatives. If the Bottzmann model analysis unit 5 has multiple simulators for obtaining price fluctuations in different 
5 ways in addition to the Bottznnann models, as will be explained below, the probability density editor 19 collects the 
probability densities from each simulator 

[0111] The foregoing is the explanation for the outline of the nnajor components of the Boltzmann analysis unit 5. 

The detailed operation flow of the Boltzmann model analysis unit 5 will be explained later again. 

[01 12] The output unit 6 presents the calculation results obtained from the system, such as the price distribution, the 

10 risk distribution, and the integrated risk index for the selected financial products. The output unit 6 is any known output 
unit as long as it can output the operation results with some forms. For example, the output unit 6 can be a printer as 
a hard copy, a display monitor as an image, or any communication means as an extemal data file. The output unit 6 
can also output intemnediate r^utts from the Bottzmann model analysts unit 5, such as the price fluctuation simulation 
results or the probability density distribution at the given trial. The output unit 6 may physically includes multiple output 

'5 means. 

[01 13] The total cross-sec^k>n/stochastk: process setup unit 7 supplies a set of sampling time intervals through setting 
ftuctuatk)n probabilities (or frequencies) per unit time at the sampling unit 14, as has been explained in conjunction 
with the sannpling unit 14 of the Boltzmann analysis unit 5. 

[01 14] TTie total cross-section/stochastic process setup unit 7 acquires the frequendes of price change and the price 
20 change rates for the selected financial product or its derivatives from the maricet database 1 1 , and supplies total cross 
sections, whk:h are defined as the ratios of the frequency of price change to the price change rate, to the Boltzmann 
equation. 

[01 15] The total cross-section of the Boltzmann equation conresponds to the frequency of price change for the finan- 
cial product or Its derivatives, whk^ will be described in more detail below. If a sampling time grid is used instead of 
25 the frequency of price change, as in the conventional syst^, then the total cross-section/stochastic process setup 
unit 7 can be omitted. 

[0116] The velocity distribution/direction distribution setup unit 8 supplies the distributions of the price change rate 
and the price change direction to the simulator 15, as has been explained. 

[01 17] The velocity distribution/direction distribution setup unit 8 acquires the past records for the financial product 
30 from the market database 1 1 , and obtains the distrft)utions of the price change rate and the price change direction from 
the past data. The distrbutions are then supplied to the simulator 15. 

[0118] Preferably, the velocity dtstrftnition/direction distribution setup unit 8 has a numerical analysis function. With 
this functk>n, the distribution setup unit 8 is capable of estimating the distribution of price change rate from the past 
data using a sigmoid function and its approximate function. Alternatively, the distribution setup unit 8 can determine a 

35 set of parameters for the sigmoid function for the distribution of the price change rate using the price change rate of 
the previous day. The distritnition s^up unit 8 also estinnates the distribution of the price change direction taking the 
conelation between the probability of price-up and the probability of price^own into account, or attematively, taking 
the correlation between the price change rate and the price change direction into account. The distribution input unit 
8 is also capable of generating a probability dlstn*bution conresponding to the price, if the distributions of the price 

^ change rate and of the price change directk>n are dependent on the price. 

[01 19] The random number generator 9 produces a random number used in the Boltzmann model analysis unit 5 to 
simulate a price fluctuation using the Monte Cario method. The generated random number is supplied to the price 
fluctuatk)n simulation unit 15, as has been described at)ove. Applk:ation of this random number will be explained below. 
[0120] The velocity distribution/direction distribution setup unit 8 and the random numtser generator 9 are essential 

^ elements for the Bottzmann model. 

[0121] The \^B-evaluation unit 10 calculates a risk or a risk distribution from the price distribution for a selected 
financial product or its derivatives. 

[0122] Several devices for cak^lating a risk or a risk distribution for a financial product or its derivatives from the 
price distribution for the financial product have been conventionally realized as VaR-evaluation device. That is, the 
50 provability density computed by the Bottzmann model according to the present invention can be transfened to any 
conventional VaR-evaluation unit 10, whbh then provides a price distribution, a risk distribution, and a risk value inte- 
grated as an index for the outputs. 

[0123] The VaR-evaluation xxivX 10 is not essential to the present inventk>n because any conventional devices may 
be used as the VaR-evaluation unit 10. 
55 [0124] Thenrtar1cetdatat>ase11 stores infonmation about ftnandal products and their derivatives. In the specification, 
the tenrn "database" includes data itself systematically stored in the database, data search engines, and hardware- 
capable of storing the data. 

[0125] Although the mark^ database 11 is essential to the present invention, any existing extemal database may 



13 



EP 1 178 416 A1 



be used. To this end, a system using an existing database together with the other essential elements of the present 
invention falls into the scope of the present invention. 

[0126] The foregoing is the configuration of the price and risk evaluation system 1 of an embodiment of the present 
invention. Preferably, the elennents conducting data processing may be included in the CPU of a computer that activates 
5 Installed programs and control^ the respective tasks. In this case, different processing means may be Included In the 
same hardware with parallel ways. The input units anrK>ng the above-described elements may be an ordinary keyboard 
or pointing device. If data Is acquired from other data files via data communication, the data communication means 
itself becomes the input unit. 

[0127] Hereafter, the theoretical backgrounds are described to clarify the reason why the Boltzmann model intro- 
10 duced into the present Invention to conduct risk analysis for financial products can con'ectly estimate the price change 
probability distributions for the financial products with higher freedom than the conventk>nai financial models. 
[0128] Basically, the present Invention needs to input several parameters, for the Boltzmann equation, which include 
an Initial price, a distribution of price change rate, a distribution of price change dlrec^on, and time domain concemed 
for each financia) product With these parameters, the Boltzmann equation is solved using the Monte Cario method, 
IS and the price and risk distributions obtained at a spedfc time domain as solutions of the equation. 

[0129] Rrst, the tinne-dependent behaviors of prk:e for a selected financial product or its derivatives are described 
by a Boltzmann equation. The present invention applies a neutron transport Boltzmann equation, whch is generally 
used to design a nuclear reactor as the established methodology in the nuclear industry. 

[0130] Neutron transport Boltzmann equation Is an equation for describing a macroscopic behavior of neutrons. A 
20 model for explaining a phenomenon based on a Boltzmann equation is calted the Boltzniann model. 

[0131] In a neutron-relating Boltzmann model, the position of a neutron is defined by the seven-dimensional vectors 
r, vQ, and t. Here, r = (x, y, z) denotes a vector in the position space of the Cartesian coordinates, vSi - (vSl^, v£2y, 
vQ^) is vector in the vekx:ity space, and t denotes time. A set consisting of the seven-dimensional vectors is called a 
phase space. 

25 [0132] A general formula of the neutron transport Boltzmann equation is as follows. 



[0133] In the equation (6), the important quantities are Z^Cr, v) and ZJjr, v*, v, £2), whk^ are a mk^-oscopic total 
cross-section and a double differential cros&-section, respectively. These quantities indk:ate probabilities of neutron 
collisions and scattering per unit length. 

[0134] MIcroscopx: cross^section is a product of an atomk: number density (unit is cnr^ and a microscopic cross- 
^ section in the reactor physk^s. The mk:roscopc cross-section is detennlned by the nudkie (for example, uranium, 
oxygen, hydrogen, and so on) existing in a nuclear reactor. On the other hand, the microscopic cross-section is an 
effective cross-section of a nucleus (unit is cm^) giving the collision probability between one nucleus and one neutron. 
The name "cross-section" s derived from the unit that expresses an area of nucleus. 

[0135] Since a concept similar to the atcmuc number density has not been found in the financial engineering field, 
^ the mk^roscopc cross-section and the microscopk: cross-section can not be distinguished. Thereupon, if the neutron 
transport Boltznnann equation is applied to finance, the microscopic cross-sectk>n and the mk^roscopic cross-section 
are consolidated in a single conc^t for cross-section. In physcs, the double differential cross-section corresponds to 
the velocity and the angle distributions of a neutron emOted from the nuclear reactions. 

[0136] It is impossible to solve an arbitrary Boltznnann equation analytically with a deterministk: way, because the 
so cross-section changes connprtcatedly in the phase space. On the contrary, it is known that the Monte Cario method is 
capable of solving the Boltzmann equation without any approximation. 

[0137] To explain the means for sohring the problem, equation (6) is simplified, and a one-dimensional and homo- 
geneous problem with no internal neutron source will be discussed. In the neutron transport theory under a one-di- 
mensional problem, a homogeneous cfistribution is assumed the y and z directions. Accordingly, the vector £2 is a th ree- 
S5 dimensional vector In the one-dimension problem. In order to discuss the one-dimensional and homogeneous problem 
applied to the finartdal fiekl, the Brownian nrotion of a price for a single stock will be described hereafter. The direction 
is denoted by \i indicating "rise" and "fair only, and the other directions are not defined. 
[0138] Equation 6 is written as equation (7). 



30 




Vdt c> . , . , , 

-Jdv'dQ'2,(r,v',Q'-»- v,Q)0(r,v',Q';t)-S(r,v,Q;t) 



(6) 
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vdt ax ^ 

5 (7) 

-Jdv'diuX(v',iu' v,^)0(x,v',Ai';t) 

[01 39] If the velocity v of equation (7) can be integrated Into the representative velocity u, and if the angle distribution 
10 Is isotropic, then equation (7) can t>e deduced to equation (8). 

9»(x:t) _p 9'»(x:9 
*3t I' 

15 

[0140] the diffusion coefTKHent becomes 



20 



here, f is the degree of freedom of the system. In this one-dimensional and homogeneous problem without an internal 
neutron source, f takes a value of '1 ' ideally, because the other directions are undefined. 

[0141] The flux expression is very convenient for a neutron transport problem. The flux expression gives many ad- 
25 . vantages to the Monte Carlo simulation. Neutron transport Monte Carlo simulation is characterized by many effective 
variarK^-reduction techniques. These techniques can be introduced using flux expression. However, describing a fi- 
nancial Monte Carlo with the flux expression likely to cause confusion, and therefore, the conventional density ex- 
pressk)n will be used temporarily. 

[0142] A neutron density function p (x, v, p.; t) is given by the solution of Boltzmann equation (1 0). 



30 
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3pOc,v^ = vn ^J^h:^ + 2:t(v)vp(x,v,H; t) 

dt dx 

- idv'dp'£3(v',|i' ^ v,nKp(xy .^i'; t) (10) 



[0143] Equation 8, whk^ is a neutron diffusk>n equation, Is also rewritten as 

^ = D'^ (11) 

[0144] The density functton P (x; t) Is an integration of p (x, v, p; t) at velocity v and angle p. The diffusion constant 
in the density expression fc>ecomes 



45 



so 



D = jj^ (12) 

here, is a frequency of collisk>n between a neutron and a medium, and is expressed as 

X, = ul,(u). (13) 



[0145] If the corrcept of the drffusion constant D' and the concept of the volatility are the sanrte, then equation 12 
stands, and the volatility o becomes 
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(14) 

[0146] Equations 13 and 14 show that the neutron velocfty and the total cross^ection are connected to volatility. A 
typical volatility a takes a value of square root of 0.1 per year [year^]. If the price changes once a day, the collision 
frequency becomes 365/year. Here, u is evaluated as being 0.0117/day, which is substantially equal to the average 
of the stock-price change rate per day. 

[0147] A numerical example of the total cross section is indicated by equation (15). 

o 

2 X 0.0117[1/dayl x 365[daysyyear] 
« 0.1[1/year] 

= 85.4 (15) 



10 



[0148] The total cross-section (£J is in inverse proportbn to square of volatility, as shown in equation (15). This 
20 relationship guarantees the equivalence between the volatility in the financial technology and the total cross sec^on 
in the Boltzmann model. 

[0149] Based on this relationship, the Boltzntann equatk>n can provide a basis for determining a pnce or a risk for a 
financial product or its derivatives by relating the variables in the neutron transport Boltzmann equation, such as position 
X, velocity v, angle \i, and time t, with the price x, the price change rate v per unit time, the price change direction \i, 

25 and the transient time t for a financial product or its derivatives, respectively. 

[0150] No concept relating to the concept of neutron double differential cross section has been defined yet In the 
financial technology. However, in order to apply the Boltzmann model to the financial technology, the double differential 
cross section should be defined. In this embodinnent, an evaluatk)n technique for a neutron cross sectk>n is applied 
directly to determine the double differential cross section of a stock price. 

30 [0151] The cross section can be evaluated from the e)q;>erimental data and theoretical computation due to nuclear 
physics for the neutron transport problem. At present, there is no established theory for evaluating the cross section 
for a stock price, and accordingly, the double differential cross section is detenmined from the market data, that is, for 
example, stock prices announced on newspapers, internets, and so on. 

[0152] In the evaluation, the distribution of the price change rate v per unit time is estimated from the past records 
35 for stock prices using a sigmoid function and the approximation form. 

[0153] In the Boltzmann model, the velocity (fistrft>ution term of the double differential cross section must be deter- 
mined by defining a sigmoid function of the price change rate v using the price change rate v* corresponding to the 
daily retum for the previous day. 

[01 54] Furthemnore, in solving the Boltzmann equation, the direction distribution term of the double differential cross 
40 section in Boltzmann equation can be detemrined as the directk>n change probability for a financial product or its 
derivatives, from the past data for the financial product. 

[0155] If there is a correlatk>n between the price-up (or upspring) probability and the price-down (or downfall)-prob- 
abilrty along the price change direction for the financial product, this con-elation can be taken into accountf or determining 
the direction distributk>n temn of the double differential cross section in the Boltzmann equation. Taking the correlation 

45 into account allows a probability drift to be evaluated in d^'l. 

[0156] There are two types of price distributions, nannety, a symnnetrical distribution with respect to the expected 
value, and an asymmetric distrtoution with respect to the expected value. A symmetrical distribution can be obtained 
by separating variables into the velocity distribution term and the direction distribution temn of the double differential 
cross section. An asynnmetrical distn^utran can be obtained by taking the correlation between the distribution of the 

50 price change rate Emd the probability of the directk>n change into consideration. 

[01 57] Under the Boltzmann model, there are two types of problenrts, namely, a homogeneous problem, in which the 
cross sectk)n is Independent of a price, and a h^erogeneous problem, in which the cross-section is dependent on the 
price. The homogeneous problem treats the distrft^utions of the price change rate and the prk:e change as constant 
concerning the price x for the financial product. On the contrary, the heterogeneous problem treats the distributions of 

55 the price change rate and the price change direction as variable in response to the price x for the financial product. 
[0158] In this nnanner, the Boltzmann model is capable of comprehensively describing both cases in which the vol- 
atility is dependent and independent of the price, by applying the Boltznrann equation to a homogeneous problem and 
a heterogeneous problem. In contrast, it was difficult for the conventional technique to treat such a realistic situation 
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with the volatilrty dependent on the price. 

[0159]. The Boltzmann model is also capable of treating both a linear problenrt, in whbh the cross section is inde- 
pendent on the probability density or the flux for the financial product, and a nonlinear problem, in which the cross 
section is dependent on the price with a consistent manner. 

5 [0160] Again, the Boltzmann model can comprehensively describe both cases, in which the volatility is dependent 
on the price as non-linear problem and independent of the price as linear problem in a consistent manner 
[01 61 ] Another advantage of the Boltzmann model ts a fact that a product (or a flux) of the probability density function 
and the price change rate per unit time for the financial product is introduced for applying a variance reduction nnethod 
in order to improve the computation efficiency. 

10 [0162] Introducing the concept of flux makes it possible to simulate a price fluctuation for a financial product at 
art)itrary time. To be more precise, the profc>ability of the price for the financial product is evaluated from the estimated 
track-length of the flux for the financial product, which is just like the track-length estimator in the neutron diffusion 
equation. This arrangement is capable of reducing the variance in the Boltzmann analysis using a Monte Cario method. 
[0163] The Boltzmann model also introduces the concept of point detector used in the neutron transport calculation, 

IS in addition to the concept of the flux, into evaluation of price fluctuation for the financial product or the derivatives. This 
arrangennent can also reduce the variance of the Boltzmann model analysis using a Monte Cario method. 
[0164] For example, using all of or a part of the events of price fluctuatk>n for the financial product to evaluate the 
probability at an infinitesimal time band and an infinrtesinnal price band, the variance can be reduced. No events of 
price change are likely to occur or no fluxes can pass through in such an infinitesimal band during the random sampling, 

^ or no fluxes can pass. 

[0165] Another method for reducing the variarK^ is to cak:ulate an adjoint probability density and an adjoint flux of 
the price change for the financial product or the derivatives, and weight the sampling in proportion to the adjoint prob- 
ability density and the adjoint flux. This method is ato based on the introduction of the concept of the flux. 
[0166] The Boltzmann model is applicable to multiple financial products or the derivatives that requires consideration 
25 of the connelation among these products. 

[0167] When evaluating a portfolio that isftnarKsal derivative product consisting of a combination for multiple financial 
products, the conrelation among the financial products or their derivatives are taken into account, the Boltzmann model 
can be adopted to the conventbnal evaluation system. 

[0168] In thefiekl of financial engineering, the tto's theory is well-known. According to the Ito's theorem, if a financial 
30 product obeys the Ito's process, then the derivative product also ot>eys the Ito's theorem. The Ito's theorem is applicable 
to the Boltzmann model as well. Accordingly, when evaluating a price or a risk for a financial derivative product in the 
Boltzmann model, the price distrft>ution of the derivative product can be propagated by the original underlying assets, 
based on the Ito's theorem. Thus, technique of the present invention can be applied to conventional price and risk 
evaluation systems. As a result, the Boltzmann model can be employed in the conventional systems. 
35 [0169] Finally, the Boltzmann model of the present invention carries out simulations for price fluctuation, and accu- 
mulates the probability distribution of individual simutatbn to obtain the prk^ distribution and the risk distribution. Ac- 
cordingly, simulations of price fluctuation can be executed in the parallel manner by the price fluctuation simulator 15 
and the probability density computation unit 18 to improve the operation speed. 

[0170] The foregoing is the theoretical background of introducing the Boltzmann model for risk analysis of the financial 
40 product. Subsequently, the price distrituition and the risk distribution for the financial product by the Boltzmann model 
will be explained in detail using actual examples of analysis results. 

[0171] Fig. 2 illustrates the operation fk>w of the Boltzmann model analysis unit 5. The Boltzmann model analysis 
unit 5 executes the steps A through t, as shown In Rg. 2. 

[0172] Step A is initialization of the Boltzmann equation. The initialization unit 1 2 executes this step. 
^ [0173] Step B gives initial values of Boltzmann equation. Solving the Boltzmann equation is equivalent to obtaining 

the Green's function of the Boltzmann equation having a source term. The Green's function r^resents virtual particle 

diffusion distribution in the phase space from a poirrt source. By setting the initial values Xq, Vq, and for the price x, 

the price change rate v, and the price change direction \i, the source temn of the Boltzmann equation is given. This 

step is canied out by the initial value setting unit 13. 
so [0174] Step C detemfiines a sannpling method. In this example, it is assumed that the price changes once a day. 

Sampling in accordance witii the frequency of price change based on the total cross-section will be described later. 

The sampling unit 14 executes step C. 

[0175] Step D simulates price fluctuation using the Monte Cario method. The integral 

55 
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^ia =J[;,' dv vexp 



V 

T 



(16) 



integrates formula v • exp[-v/r|, which form has similarity with an empirical fonm for neutron kinetic energy emitted in 
fission reaction with the velocity temn in the differential cross-section. In other words, price fluctuation is simulated by 
giving a velocity distribution to the Bottzmann equation, and by calculating v^^ satisfying the probability associated 
10 with random numt)ers having been generated. The price flucmjation simulator 15 carries out step D. 

[0176] In the example shown in Fig. 2. the price change direction is simulated In such a manner that the price-up 
direction equals 1 (^^4.1=1 ) and the price^own direction equals -1 0H4.i='1 ) simply from the value of the random 
number. 

[0177] Step E carries out integration of the Green's function to obtain the probability density Pm. The probability 
'5 density computation unit 1 6 carri^ out step E. 

[01 78] Step F (tetemnines whether or not one trial has been completed. The one-trial completion detector 1 7 executes 
step R If one-trial has not been completed yet, the proc^ returns to step C. 

[0179] Steps G and H d^ermine whether or not all trials have been completed, tf all trials have not been completed 
yet, the process returns to step B. The all-trial completion detector 1 8 carries out steps G and H. 

20 [01 80] Step t edits the calculated probability density, which is carried out by the probability density editor 1 9. 

pi 81] In the example shown in Rg. 2, sampling is conducted on the assumption that a price change occurs once a 
day. Rg. 3 illustrates another example. In which the sampling intenml is set in response to the frequency of price change. 
[0182] The process shown in Rg. 3 is almost the same as that shown in Rg. 2, except for steps C and E*. In step 
C, the sampling interval set in accordance with the frequency of price change. The microscopic total cross-section in 

25 the Boltzmann equation means the inverse of the mean free path (which is the average distance from one collision to 
the next collision) of the neutron. A collision frequency (which is the collision probability per unit time) is the product 
of the total cross-section and the velocity v. By applying this to the change of a stock price, simulation of price fluctuation 
can be carried out without setting a time grid, as long as the stochastk: process of the change in the stock price and 
the total cross-section are known. 

30 [01 83] In the conventional technque, using At as a certain change in time, the price after At is simulated by a nomrtai 
random numt)er having a standard dev^on of cW At. For this reason, the time grid At must be set as an indispensable 
step for conducting an accurate simulation. On the contrary, with the price and risk evaluation system of the present 
invention, setting of the time grid can be omitted in step C*. 

[0184] In step C* shown in Fig. 3, the index distribution used in the simulation of price fluctuation is expressed as 

35 

^_ -Nl-^i) 



40 [0185] In step C, the stochastic process and the total cross-section are supplied from the all cross-se<^on/stochastic 
process input unit 7. 

[0186] In step D of Rg. 3, price fluctuatbn is simulated in accordance with the sampling method mentioned above. 
The price fluctuation simulation itself is substantially the same as the simulation of step D shown in Fig. 2. 
[0187] The difference is that the sampling interval changes in response to the frequency of price change in the 
45 process shown in Rg. 3. The sampOng interval is adjusted by ctetenmining whether the next sampling position reskies 
within the obsenfatk^n area (Am, Bm, Cm, Dm) after every price sirmjlatk>n. 

[0188] Rgs. 4 and 5 scherrmtkrally show the simulatk>n within a predetermined obsenratk>n area corresponding to 
the processes of Rgures 2 artd 3, respectively. Espedalty, it is dearly shown in Rg. 5 that the price f luc^uatbn can be 
simulated at a spedfk: time interval depending on the frequency of price fluctuation, irrespective of the time grid. 
50 [0189] The simulation method of price fluctuation based on the Bottzmann model according to the present invention 
can also carry out simulatk>n based on the conventional diffusion model if equation 1 8 is used as a function describing 
time T, directk>n M, and velocity V. 
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|pi90] The results of two examples will be shown below. 
<Example 1 > r = 0.05. = 0.1 1 , and T = 0.25 year. 

[01 91] The lower and upper limits Am and Bm of x (price) shown in the flowchart of Rg. 3 define a price band, which 
is determined by dividing the range of -3a<x<3o (Sx = 0.1a) by 60. The lower limit Cm and the upper limit Dm of time 
t define a time band, where Cm equals 0.25 [years] and Dm equals CrrH-( 1/365) [years]. An evaluation quantity a)| is 1 . 
[0192] Fig. 6 shows the evaluation result of example 1 using the solid line 22 in comparison with the theoretical 
distribution (i.e., the logarithmic normal distribution) indicated by the dashed line 21 . As is dear from fig. 6, the simu- 
lation result of the present invention indicated by the solid line 22 is almost coincident with the theoretical distribution 21 . 

<Example 2> = 0.1 , r = 0.05 and r = 0. 

[01 93] The lower and upper limits Am and Bm of x shown in the flowchart of Fig. 3 are set to Am = and Bm = +»>. 
The time band defined by Cm and Dm are from 0 to 365 [day] with 5| = 1 [day]. The evaluation quantity o)| Is x. In Fig. 
7, the dashed line 23 indk^ates the theoretical distribution 23 under a drift, and the long dashed line 24 indicates the 
theor^ical distribution 24 without a drift. The simulation results 25 and 26 obtained in example 2 substantially reproduce 
the theoretical distributions with and without a drift. 

[0194] The velocity distribution and the direction distribution used in the above-explained simulation are the proba- 
bility distribution same as the standard Gaussian distr^ution. Consequently, the probability density obtained by the 
simulation becomes equivalent to that of the diffusion model. 

[0195] In order to realize the Boltzmann model, the velocity distribution and the direction distribution must be eval- 
uated. 

[01 96] An example of how to evaluate the velocity distribution will be now explained. To apply the geometric Brownian 
motions, the natural logarithm of stock price is defined as x. This corresponds to the position x of a neutron. 
[01 97] From stock prices of about sixty Japanese electrk: m^inery makers over three years, the vekx:fty distribution 
and the direction distribution are evaluated. To define x, the natural togarithm of the closing price (or the last price) of 
each day is input to x. An irtcident vekx:ity v* is defined as the absolute value of the difference between natural logarithm 
of the dosing price of the cunBnt day and natural k^garithm of the dosing prk^e of the previous day. A current velodty 
V is defined as the absolute value of the difference between the natural bgarfthm of the closing price of the current 
day and the natural togarfthm of the dosing price of the next day. 

[0198] The inddence direction (i* is represented by the negative or positive sign of v*, and the current direction \i is 
represented by the negative or positive sign of v. Since a change per day is observed, the deterministk: drift temn (e. 
g., a non-risky interest rate) is omitted. If any drifts are found In the simulatk>n using the actual data mentioned above, 
it is a purely stochastic drift. 

[0199] In order to obtain the velocity distribution of the stock price, the spectrum shown in Rg. 8 is required. The 
spectrum is the integral of x, ^. t of the density p (x, v, \i; t), and is expressed by equation (19). 



S(v) = Jdtd^dx * p(x,v,^; t) 



(19) 
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[0200] In Rg. 8, the darkened drde (•) indicates the total spectrum 27 expressed by equation (19). The spectra of 
the negative (or price-down) direction and the positive (or price-up) direction are indicated by stars (•) and white squares 
(□) 29, respectively. The price-down spectrum S_ (v) and the price-up spectrum (v) are expressed as 

S_(v) =. fdtdx • p(x,v,-l;t) 

; (20) 
S,(v)-Jdtdxp(x,v,l;t) 

[0201] Although the two spectra are slightly different in size, the shapes are the same. This resuK shows that the 
velocity distribution and the direction distribution are Independent of each other. These distributions are described by 
the sigmoid function expressed by equation 21 artd its approximation using the Maxwell distribution as a typical exam- 
ple. 

^iL^^f^r^ (21) 

^ gexp[Tfv] + ;-g 

here, v, k, g, and yare arbitrary real numbers. 

[0202] These spectra are approximated as indicated by the steep slop 30 and the gentle slop 31 . In other words, the 
two spectra correspond to the two components* namely, the steep slope component and the gentle slope component. 
The curve 32 indicates the Giaussian distribution. The Gaussian distribution almost reproduces the steep slope com- 
ponent, but it evaluates the gentle slop component excessively small. 

[0203] Rgure 9 exhibits the dependency of the spectra on the incident velocity. The darkened marks 33, the cross 
marics 34, and the white martcs 35 represent the velocity distributk>ns with the incident velocities of about 1%, about 
2%, and about 3%, respectively. These distributions are normalized to 1 .0 with integration. 

[0204] Rgures 1 0 and 1 1 Illustrate the double drfferential cross-section £ (V, m.'-v,^) with respect to direction \jl From 
Figs. 1 0 and 1 1 , it is apparent that the shapes of the spectra are the same. This fact rndk^tes that the double drfferential 
cross-section is given by the product of the velocity distribution 3 CV->v) and the direction distribution p '(^'^). This 
can be expressed by equation 22. 



2(V', A*' ^ V,Ai) OC S(V' -V) PifJL' - fM) (22) 

[0205] Here is an example of determining the distribution of the price change by a sigmoid function, using the past 
records for price change rate as a pararrteter. As Illustrated in Rgs. 9 through 11 , the spectra shift to higher values as 
the incident velocity v increases. Therefore, the concept of temperature can be introduced, and the distributions shown 
In Rgs. 9 through 11 can be described by the Maxwell distr&ution. 

[0206] Since the distributions shown in Rgs. 9 through 11 suggest an exponential distribution, the evaporation spec- 
trum used to describe the neutron emitted from a nuclear reaction, whbh is expressed by equation (23), will be used. 



f(v) OC vexp 



T 



(23) 



[0207] This is a modification of the Maxwell dt^ributbn. The natural logarithm of equation 23 becomes 



lntf(v)) - ln(v) = - ^ + Const. 



(24) 



[0208] Rg. 12 illustrates the relationship of equation (24). The inverse of the slope corresponds to temperature T 
Rg. 12 gives an experimental equation of differential cross-section expressed by equation (25), and Fig. 13 illustrates 
the relationship between the velocity v* and temperature T. 
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3.5533V 0.1023V '+ 0.0044 



(25) 



[0209] Next, an example of estimating the probability of the change direction for a financial product or the derivatives 
from the past records of that financial product will be explained. 

[0210] In this problem, the direction takes values of only 1 and -1 . The value "1 " denotes increase in price, and "-1 " 
means decrease in price. In the financial engineering, the direction distribution is given by equation (26). 



continuously price- up 
change f r oinpr i ce- dovn t o pr i ce- up 
change f r ompr i ce- up t o pr i ce- dom 
continuously price-dow 

(26) 

[0211] Fig. 14 illustrates the averages of the continued price-up probability (1 -> 1 ; t) and the continued price- 
down probability (-1 -1 ; t) of every five days. The daricened squares (■) 36 represent the events transient from 
price-up to price-up, the probability of which is expressed by ^< (1 1 ; t ), and white squares (D) 37 represents the 
events transient from price-down to price-down, the probability of which is expressed by (-1 - 1 ; t). The bold 
horizontal solid Iine38 and the dashed line 39 are the time averages of these two probabilities. Other two probabilities 
are expressed by equation (27). 

p(i-l:Ut)-i-p(i-itt) 

-l;t)-l-p(-l- -l;t) 

[0212] Fig. 14 exhibits the conretation between the probability of price-up and the probability of price-down with 
respect to the probability of the change direction for a financial product or its derivatives. Rg. 14 clearty shows that 
probability of price-up (denoted by ■36) and the probability of price-down (denoted by □ 37) change in opposite 
directions as time passes. This fact indicates a negative correlation. 

[0213] By reflecting the conretation shown in Rg. 14 into the angle distribution of the double drfferontial cross-section, 
a more accurate evaluation of the stochastic drift can be achieved. 

[0214] Next, connparlson will be made betwe^ the evaluation results with the price and risk evaluation system of 
the present invention and the evaluation results with the conventional technique. 

[0215] Fig. 15 shows the evaluation results with the Boltzmann model that uses the price change rate distribution 
and the price change direction distrbution. The solid lines 40 represent the results from the Boltzmann model, which 
effectively reproduce the jumps (big changes) in price appearing in the thicic line 41 that indicates the actual record. 
[0216] Fig. 16 is a detailed view of Rg. 15. The track depk:ted by the darkened squares 42 is the real record. The 
real record exhibits several big changes (jumps) 43 of about 10% per day. 

[0217] The white squares (□) 44 and the white triangles (A) 45 are obtained from the simulation with the Boltzmann 
model. These symbols exhibit jumps 46 that are sbnilar to the jumps 43 in the real record 42. The ability of simulating 
price jumps is the significant feature for the price and risk evaluation system of the present invention. 
[0218] I n contrast, the conventional diffusion model (denoted by cross marics 47 in Fig. 1 6) is incapable of reproducing 
the abmpt jumps, and it simulates price fluctuation only in the continuous manner. 

[0219] Fig. 1 7 illustrates the simuladon result obtained by the Boltzmann model in comparison with logarithmic dis- 
tribution of a stock price of after 200 days. In this simulation, the deterministb drift temn (for example, a non-risky 
interest rate) is not taken into account. In Rg. 17, the solid curve 48 indicates the conventional diffusion model (with 
the expectation value of 0 and the standard deviatk>n of 0.29). The darkened squares 49 indicate the Boltzmann model, 
whk:h is slightly drifted due to purely the stochastic process. The dashed line 50 indk:ates the con-ected Gaussian 
distribution as a result of correction of the drift tenm. The corrected Gaussian distribution 50 covers a range of at>out 



P(i-i;t) 

|^(-l-l;t) 

P(i--i;t) 

p(-l--l;t) 
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±2.5a of the Bottzmann model; however, the other portions were underestimated. 

[0220] The simulation result shown in Fig. 17 exhibits the significant feature of the Boltznnann nrK>del well. That is, 
the Boftzmann model used in the present invention is capable of evaluating the stochastic drift and the probability of 
big jumps in price, which can not be systematically reproduced by the conventional technique. 

5 [0221] Next, another application will be explained with a case in which a conflation between the velocity v (conne- 
sponding to the price change rate) and the angle \l (corresponding to the price change direction) exists. In the above- 
described examples, the velocity v and the angle \l are supposed to be independent of each other with respect to the 
double differential cross-section, and therefore, separation form of variables was applied, as expressed in equation (22). 
[0222] However, as the significant feature, the Bottzmann model is capable of taking the correlation between the 

10 velocity v and the angle \i into account by introducing a function that is not subjected to separation of variables. The 
above-described examples without considering the correlation only evaluate synrimetricai distributions, which are sym- 
metric with respect to the mean value. On the contrary, in this example, an asynnmetric distribution can be evaluated 
by taking the correlation between the probability of the prk:e change rate and the probability of the price change direc- 
tion, and therefore, by introducing a furYction of non-separation of variable. 

IS [0223] Next, a heterogeneous problem will be explained. In equation (7), the cross-se(^on Z is constant with respect 
to the prk:e x. This is the same thing as the conventional financial engineering in whk:h the volatility is constant with 
respect to the price. When the conventional financial engineering treats a heterogeneous problem containing an in- 
constant cross-section £, the price-dependency of the volatility had to be corrected by a technique of volatility smile 
or other techniques. However, these techniques greatiy rely on past experiences and know-how. In contrast, the Boltz- 

20 mann nrKxlel applies equation (6) expressing a heterogeneous problem , which allows the price-dependency of volati lity 
to be systematically considered. 

[0224] If the heterogeneity varies due to the change in the price distribution, a non-linear Boltzmann equation will 
be introduced. Since the conventional financial engineering is not able to theoretbally treat heterogeneous problems, 
these problems are often treated as stochastic volatilities in the conventional technique. However, stochastic volatility 
25 itself greatiy relies upon know-how or experiences, and is lack of objectiveness. On the contrary, the Boltzmann model 
. used in the present invention can treat stochaslrc volatility accurately using a non-linear Boltzmann equation (28). 

-£$5^ = ^ ^*?^ + E,W(x.v,^;t).vW(x.v.^:t) 

vdt dx 

30 

-Jdv•d^•2:^((|►(x,v,^i; t)v*.ji' ^ v,^)<|>(x,v',M.'; t) (28) 



[0225] In equation (28), flux <^ is contained in the cross-sectk>n, and the cross-section changes along with the change 
in flux. This technique s actually used in computation of buming of nuclear fuel in the neutron-relating field. By applying' 
this technique to the financial engineering, a systematb evaluation method for stochastic volatility is realized. 
[0226] Next, the evaluation at an arbitrary point of time within an observation area will be explained. The density 
expression in equation (10) can not evaluate the probability density of an arbitrary hour between the 199th day and 
the 200th day because the event of price change can not be detected. 

[0227] However, ^ the flux expressions of equations (6) and (7) are used, flux <^ (x, v, ^; t) is obtained at an arbitrary 
time t inrespective of presence or absent of price change. Since flux ^ (x, v, \i; t) can describe an arbitrary time t, and 
therefore, the probability density P at an arbitrary hour between the 1 99th day and the 200th day can be obtained 
correctiy using equation 29. 

P = /m=1*^^ (29) 

P22SI This is a track-length estimator. 

[0229] Next, a method, for nnaking evaluation at a point that does not allow effective sampling, will be explained. 
Although the estimation of the track length allows evaluation at an arbitrary time, it is impossible for the conventional 
technique to effectively collect sampling data at an intinitesinnal price band and an infinitesimal time band even if in- 
tensive sampling is carried out. This nneans that evaluation can not be made from the estimation of track length. 
pi230] In the present invention, the concept of point detector for evaluating a neutron at an arbitrary point in the 
phase space is applied to the events for financial products in order to allow evaluation at an infinitesimal interval. With 
the concept of point detector, a neutron that reaches point C is evaluated by computing the probability that the neutron 
starting from point A collides at point B and is scattered toward point C. The probability that the neutron passes point 
C is estimated from scattering infomnation that the neutron does not pass point C. In the simplest example, with the 
distance r between points B and C, the neutron decays by exp (-£( - r), and the solid angle changes in accordance with 
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the distance. In this case, the neutron that changes at point B and reaches point C is estiniated accurately to the 
considerable extent by connection of Because it is known that the probability of the scattering angle at point B is 
the differential cross-section (v',^' -> v,^), the probability of the neutron that changes at point B during the sampling 
and does not reach point C can be estimated. 

[0231] The concept of point detector is introduced into evaluation for a financial product, white using all of or a part 
of the evertts of the price change for the financial product or its derivatives. This arrangement allows the price distribution 
or the risk dstribution for the financial product in an infinitesimal observation area (or a target area). 
[0232] In reality, no events of price change are likely to occur during the random sampling, and no flux pass through 
in an infinitesinnal price band and at an inf inrtesinr^ time. In spite of the fact that an event can not occur in an infinitesimal 
area in the target phase space, the probabifity of events in such an infinitesimal area can be evaluated by the present 
inventk)n within realists computation time by automatically checking the route of causing the events. Evaluating the 
probability in an infinitesimal area, through which flux can not pass, allows the variance to be reduced efficiently. 
PS33] Next example shows a case in whch an adjoint probability density or adjoint flux of the Boltzmann model is 
introduced. The adjoint equation of equation (7) is expressed by equation (30). 



The adjoint flux (x, v, \i\ t) is a solution of equation (30), and it represents the sensitivity of the expectation value of 
the ajoint radiation source S' (x, v, t). The expectation value is expressed by equation (31) 



[0234] The adjoint source S* (x, v, \i; t) conresponds to the price evaluation equatlori for a financial product or its 
derivatives. 

[0235] By providing weights proportional to the adjoint flux (^* (x, v, p.; t) in the phase space, the variance inevitably 
accompanying the Monte Cario method can t>e reduced when evaluating the expectation value of the financial product 
or its derivatives. 

[0236] Next, an application will be made to the case in whk:h tfiere is a correlation anoong financial products. For 
example, there is a correlation expressed by equatk>n (32) among multiple (e.g., two) financial products and their 
derivatives. 



If there is the correlatk)n of equation (32) among the financial products, the conventional technique simulates a prk» 
by generating a correlative random number in accordance with a known correlation coefficient when producing nonnal 
random numbers ^i and ^ with respect to two Ito's processes. This conventional method is capable of evaluating not 
only the price for a single financial product, but also the price of a portfolk) consisting of a combination of multiple 
financta) products. 

[0237] In contrast, the pr^ent invention realizes application to the portfolio by simultaneously setting multiple equa- 
tions (33) for multiple financial products based on the Boltzmann model. 




(30) 



-Jdv'dA*'2.(v,A4-*v',/*>*(x,v',m';t)-S*(x,v,/i;t) 



<S*{x,v.»i; t)<t>(x,v,n; t)) = JdvdndxdtS*(x.v,|t; t)(|>(x,v,fi; t) 



(31) 




(32) 



(x^.v^.ji^:t) 



+ J^(Vi)<|>i(Xi,Vi.ji^;t) 
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- JdVidji'^Z^Cv!,, Hi -> (<t»i(Xi, v!,. n!,; t) + *2(Xi , v!, , ji!, ; t)) 




- JdV2dn2i;s(v2. H2 ^ ^2) WiC^a- V2. t) + <>2(X2. V2. j4; t)) 



(33) 



[0238] By considering the connelatbn between the double differential cross-sections (v^ -> , p^) and 1^ (V2\ 
-> V2, for the third temns of the right-hand-side of the respective equations, an effect equal to or higher than the 
conventional technique using equation (32) can achieved. 

[0239] Next, the application of fto's theorenn into the Boitzmann nnodel of the present Invention will be explained. Ito's 
theorem defines that if the price S of a financial product, such as a stock, obeys the Ito's process expressed by equation 
(34), then the price F (S, t) of the derivative product moves in accordance with the stock price, and also obeys the Ito's 
process. 



[0241 ] The conventional technk^ue that does not apply th e Ito's theorem Is based on the assumption that the random 
number ^, has a nomnat distribution. On the contrary,, the Ito's theorem stands in the Boltzmann model even if the 
distribution is not Gaussian, as long as the random process of the second term is proportional to the square root of 
the infinitesimal tinrte dt. 

[0242] For this reason, the Boltzmann model applying the tto's theorem can evaluate a price distribution in^espective 
of whether or not the distribution is Gaussian. If the variance of the distribution evaluated by the Bottznnann model is 
proportional to the square root of time, the Ito's theorem can be applied to a conventional price/risk evaluation system. 
In this case, the conventional system applying the Ito's theorem will achieve the similar effect as the present invention 
by replacing the normal distribution of the random nunnt>er x with the distrbution obtained by the Boltznrtann model. 
[0243] Although a price distribution obtained by the Boltzmann model is different from a price dtstiibution using Gaus- 
sian in the strict sense, the distribution obtained by nomrtalizing the price to the standard deviation beconrte constant 
independent of time, as shown in Rg. 18. In Ftg. 18, the darkened square (■) 51 , the darkened triangle (A) 52, and 
the darkened diamond (♦) 53 show the price distribution of after 20 days, after 60 days, and after 1 00 days, respectively. 
The white square (□) 54, the white triangle (A) 55, and the white diamond (0) 56 show the price distribution of after 
120 days, after 160 days, and after 200 days, respectively. These data coincide with the normal distribution 57 within 
the range of about ±2.5 o. The standard deviation of this normal distribution is proportional to the square root of time. 
This means that the standard deviation of the prok>ability distribution obtained by the Boltzmann model is in proportion 
to the square root of time. 

[0244] The example shown in Rg. 1 8 is the simulation result of the stock-price distribution for about sixty Japanese 
electric machinery makers. If a conventional system is designed based on the Ito's theorem, it can evaluate a price 
distribution or a risk distrftiution for derivatives deriving from the stock prices of these makers. In this case, the con- 
ventional system using the Ito's theorem can replace the normal distribution with the Boltzmann model distribution in 
order to improve the predctran ability for a price distribution of the derivatives. 

[0245] Lastly, an application of the present invention to a parallel processing system will be explained. The present 
invention uses the Monte Carlo method for the numerk:al cateulations. It is widely known that the Monte Cario method 
is an advantageous technique because the processing speed can be drastically improved by parallel processing. Es- 
pecially, it ts well known that parallel processing is quite effective in application of the neutron transport Monte Carlo 
method. Since the present invention makes use of the neutron transport Monte Cario technique in the Boltzmann 
model, the cateulation speed can be effectively improved using parallel processing computers. 
[0246] Fig. 19 illusti^es an example of a parallel processing system. In Rg. 19, each operation flow A through I is 



dS = a(S, t)dt + b(S, t)M^ 



(34) 



[0240] In this case, the price F of the derivative product is expressed by equation (35). 




(35) 
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the same as that shown in Rgs. 2 and 3. Espedalty, simulation process consisting of steps B through F is divided into 
multiple parallel flow in order to allocate the trials to a plurality of CPUs. The operation speed increases depending on 
the number of the CPUs used in the parallel processing. 

[Q247] The foregoing is the preferred embodiment of the first feature of the present invention, which is realized as a 
price and risk evaluation system for a finandal product. The present invention is also realized as a computer-readable 
recording medium storing the price and risk evaluation program, which controls a computer system to carry out the 
process described above. In this case, the program Is installed in a computer system, and the price and risk evaluation 
system is realized when the program is started on the computer system. 

[Q248] The second feature of the present invention, which is embodied as a dealing system, will be described with 
reference to Rg. 20 and the subsequent drawings. First, option-pricing methodology will be explained In the frame 
wori( of the theoretk^al and realistic aspects. A price of European option for stocks and stock indices at delivery dates 
can be evaluated using integral equations (36) and (37) based on the risk-neutral probability measure (probability 
density) P (S, t) for underiying assets in the free maricet. 

Call(K,r) =. e""^ J^dSP (S,r) (S-K) (36) 



Put(K,T)==e-" dSP(S,r) (K-S) (37) 

here S is the price of the underlying assets, t is the period to the maturity, r is the non-risky interest rate (that is, the 
money rate fixed up to the maturity), and K is the exercise price. 

[0249] Equation (36) expresses the theoretk»l call option price to buy the underlying assets (i.e., call option) at the 
maturity with the exercise price K. Equatbn (37) expresses the theoretical put option price to sell the underlying assets 
(i.e., put option) at the maturity with the exercise price K. A purchaser of these options can exercise the right at the 
exercise price K irrespective of the actual price of the underlying asset at the maturity. For example, the purchaser of 
a call option can buy that option at price K, even if the underlying price (i.e., the price of the underlying assets) is higher 
than the exercise price at the maturity. The purchaser of a put option has an obligation of selling the option at the 
maturity at price K. However, this purchaser can repeatedly trade the underiying assets in response to price changes, 
and can sell the option at price K at least at cost of equation (36). 

[02^] Black-Sholes equatk>n (BS equation) is often used to evaluate an option price. If the logarithm nonnal distri- 
bution expressed by equation (38) is input to the risk neutrality probability measurement in equations (36) and (37), 
then Black-Sholes equatk>ns (39) and (40) are obtained. 



CalI(K,T) = SN(di) - Ke"" N(d2) (39) 
Put(K,T) - -SN(-d,) + Ke'" N(-d2) (40) 

here d^ and d^ are expressed by 




25 



EP1 178416A1 




[D251] The parameter a in equations (38), (39) and (40) is the price change rate (or votatillty), and is the diffusion 
10 constant of the geometric Brownian motion model, in which the underlying price diffuses with respect to the logarithm 
of the price, for the undertytng assets. 

[Q252] The Black-Sholes equation is derived on the assumption that the volatility a is constant with respect to t and 
S. Accordingly, the Black-Sholes equation assunnes a statistic market that exhibits a constant inrespective of time and 
price. 

IS [0253] However, the real market changes as time and price chartge. Rg. 20 illustrates the price change rate CI for 
the underlying assets predicted by the geometry Brownian model, in comparison with the change rate of the closing 
price (i.e., the daily earning rate) C2 for a typical stock price. Although the volatility of the two data are almost the same, 
the appearances of the price change quite differ from each other. The geometric Brownian motion model CI does not 
exhibit a big price change, whereas the actual market price stgnrficantty varies as indicated by the curve C2. This 

20 comparison result leads to the conclusk>h tttat it is diffk:ult to evaluate the option price based on the Black-Sholes 
equation, if the underlying assets is the individual stock prk». In reality, the transaction of indivkiual stock option is 
small in number. 

[0254] The stock index, that is, the corrected average of the stock prices of many issues (for example, the Nikkei 
225 Stock Average) moves more moderately than the indivkJual stock price. Accordingly, it becomes easier to evaluate 

25 the option price forthe stock based on the Black-Sholes equation, and many transactions are canied out at the present. 
However, even if taking Nikkei 225 stock average, the appearance of the daily eaming rate C3 is still different from the 
geometric Brownian motion model (curve CI in Rg. 20), as shown in Rg 21 . If curve C3 is compared with curve C2 of 
actual record shown in Rg. 20, these two curves are essentially the same, except for the size in change. To this end, 
dealers of stock index optk>ns generally use modified Black-Sholes equations to evaluate option prices. 

30 [0255] Rg. 22 illustrates an example of the correction. The implied volatility (IV) is defined as the volatility implying 
an option price actually traded in the maritet using the Black-Sholes model. The daricened squares (■) M51 shown in 
Rg. 22 represent a typical implied volatility of the closing price for a stock index put option of the Nikkei 225 Stock 
Average. 

[0256] The horizontal line C4 extending at 30% of the vertical axis is the historical volatility calculated from the motion 
35 of the option of the Nikkei 225 Stock Average. If the market completely obeys the geometric Brownian movement model 

that is the tesis of the Black-Sholes equation, the darkened squares (■) M51 should be located on the 30% line C4. 

However, in reaTity, as the exercise price separates from the underiying price, the implied volatility tends to increase. 

The tendency that the implied volatility increases from the point at whk:h the exercise price and the underiying price 

are equal, that is, with the ratio of (exercise price)/(undertying price) being 1 .00, is called a smile curve, whk:h is 
^ indcated as C5 in Rg. 22. tt is known that the smile curve has a term structure iri whk^h the curvature becomes gentle 

as the tenm (or period) increases up to the nrtaturity. 

[0257] In general, option-dealers try to grasp the volatility matrix, which bring the smile curve and the tenm structure 
of the implied volatiFrty together based on the transaction price in the market. They determine the option price by 
correcting ttie option price obtained by the Btack-Shotes equation using the volatility matrix. 
[0258] Although the volatility matrix is one of the most successful tools for evaluating an option price, it still has some 
drawt>acks. The rrtajor drawbacks of the volatility matrix are the following two: 



1 . If there is no or a few transactions, it is impossible to obtain implied volatility. 

2. The volatility matrix can not specify a typk»l transaction in the market, in whk:h the transaction price varies widely. 

50 

[0259] The drawt)ack 1 is the essential problem concerning the implied volatility, and can not be solved by the implied 
volatility. On the other hand, the drawback 2 could be solved by a filtering technique for extracting significant information 
among from the wkiely varied infonnation in order to spectfy the realistic transaction. 

[0260] Prior to applying the filtering technique, the mechanisms, why the smile and the term structure of the implied 
S5 volatility appear, must be ctarified in order to grasp the average behavior. Although the mechanisms of the smile curve 
and the temi structure have not been completely clarified yet, various researches suggest that the major reasons are 
relating to the Leptokurcity and the Fat-tail occurring in the probabilrty distributions for the actual price changes in free 
trade maricets. Leptokurcity is the phenomenon that the probability of the actual price changes obsen^ed is sharper 
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than the normal distn'bution assumed in the Black-Sholes equation in the region of small price change. Fat-tail is the 
phenomenon that the probability of the actual pnce changes widens toward the end in the region of big price changes. 
[Q261] Fig. 23 illustrates an example of these phenomena. The di^bution of the white squares (□) M52, which 
represent the actual daily earning rates, becomes sharperthan the normal distribution C6 near the center, and broadens 

5 towards the encte. Rg. 24 illustrates the probability of price change estimated from the Boltzmarin model, in comparison 
with the logarithmic normal distribution of the Black-Sholes equation. Under this price distribution, the width (that is, 
the volatility) of the probability density distribution is smaller than that of the logarithmic normal distribution C7 around 
the relative price of 1 .0, as is indicated by darkened square (■) M53. The volatility of the probability density distribution 
becomes larger than that of the logarithm^ nomial distribution C7 in the ranges of the relative price of above 2.0 and 

10 below 0.5. It is supposed that the price distribution (indicated by white squares (□) M54 approaches the normal dis- 
tribution as time elapses according to the central limit theorem. Accordingly, the peak portion and the skirt of the price 
distribution becon^e very similar to the normal distribution C8 as time passes. This is supposed to be the factor of 
emergence of the smile curve and the tenn structure. The similar discussion is made in "John C. Hull, "OPTIONS, 
FUTURES & OTHER DERIVATIVES, Fourth Edition", Prentfce-Hall International Inc., 2000, chapter 17. 

IS [0262] The Fat-Tail of the priceK:hange distribution corresponds to the big price changes that occur in the real price 
fluctuatk>ns C2 and 03 shown in Figs. 20 and 21 . There are two models taking such big changes into account, namely, 
a Jump model and a probability volatility model. The jump rrKKlel reproduces the Fat-tail independently in the stochastk; 
process that is totally different from the normal distrbution. In the probability volatility model, the standard deviation of 
the normal distribution (that is, the volatility) fluctuates wrth time. The jump model is based on the assumption of dis- 

20 continuous price changes, while the probability volatility model is essentially a non-linear problem. For this reason, 
either model is incapable of achieving the risk-neutral probability measure uniquely. Consequently, equations (36) and 
(37) of evaluating option prices can not be applied to these two models, which is the major drawback. 
[0263] In contrast, although the Boltzmann model proposed by the present invention covers the category of probability 
volatility model in a broad sense, a linear Boltzmann equation can reproduce the characteristics of the Leptokurcrty 

^ and Fat-Tail. If the angle distribution of the linear Boltzmann equation is isotropk: (that is, if (m.'->M-; t) = 1/2 in equation 
(26)), the solutton becomes risk-neutral and unique. Therefore, applying the Boltzmann model to evaluation of option 
prices allows the esserrtial trend of the volatility matrix to be property estimated. 

[0264] One of the significant features of the Boltzman n model is that the present model can treat market-dependency 
of price fluctuation. The market-dependency means that a set of big price changes occur coincidentatfy with certain 
30 time intervals. The price evaluation system that has been described above as the first feature of the present invention 
preferably recommends applying an evaporation spec^m equation (42), which is a modifk^ation of the Maxwell's 
distribution, as the price distribution f(v) in order to taking Leptokurdty into account. 

f (v) « vexp 



^ [0265] The Boltzmann model treats the correlation between the price charige rate in the underlying assets and the 
previous price change rate. The Bottzm£mn model claims the existence of a definite market-dependency between the 
daily eaming rate v* of the previous day and the daily earning rate v of the current day via temperature T as exemplified 
in Eq.(42) in case focusing on the cbsing prices. Rg. 25 illustrates a typk:al example of the market-dependency. In 
Rg. 25, the darkened squares (■) M56 represent the temperature obtained from the real records of the closing price. 

^ The curve 09 is a fitting line of the darkened squares with a quadratic function. The fitting line exhibits the fact that the 
temperature T has a quadratic tendency with respect to the daily eaming rate v* of the previous day expressed by 
equation (43). 



T(v') = To(l + CoV'+g„v'=') (43) 



[0266] The quadratk) dependency recalls a direct analogy to the instability of the stock market in a system with a 
55 positive feedback such that the specific heat increases as the temperature r^es. 

[0267] The curve 05 extending along the real records (i.e., the darkened squares (■) M51) in Fig. 22 exhibits a 
volatility smile. This volatility smile is obtained by evaluating the option pric^ of equations (36) and (37) based on the 
Boltzmann model and plotting the volatility of the Btack-Sholes equatk>n that become equal to the evaluation result. 



(42) 
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Rg. 26 illustrates the daily return rate CI 0 obtained by the Bottzmann model in the simulation process for price eval- 
^uation. The curve C 10 very well reproduces the daily earning rate (i.e., the darkened squares (■) M56) of the Nikkei 
225 Stock Avera^ in the same temn as the option transaction. 

[0268] The daily earning rate shown in Rg. 26 show^ a typical Fat-Tail. Based on this daily earning distribution 
5 obtained from the BoRzmann model, a random number ^ is generated to simulate the tracks of the underlying assets 
S using equation (44). 

f = rdt+^ (44) 

10 

[0269] The jump model descrbed above does not take the market-dependency into account, and it treats a big prk^e 
variation as abrupt and discontinuous jumps. This jump model may appear to give the simitar results as the Boltzmann 
model; however, the result is quite different from those obtained from the Bottzmann model. 

[0270] Rg. 27 illustrates the implied volatility of the Boltzmann model in comparison with that of the jump model. The 
solid line C11 and the broken line CI 2 show the result of Boltzmann model with the temis toward the maturity being 
40 days and 80 days, respectively. Another broken line C1 3 and the dashed line CI 4 show the result of the jump nwdel 
with the terms toward the maturity being 40 days and 80 days, respectively. 

[0271] The comparison result shows that the implied volatility of the jump models C1 3 and C1 4 become larger than 
that of the Boltzmann models C11 and C12, artd that the curvature of the smile curves of the jump model is snnaller 

20 than that of the Bottzmann model. This result, that is, the excessive implied volatility of the jump mode), also applies 
to the comparison between the junrtp model and the actual past records. Because in the jump model the size of a price 
change is not correlated at all, the central limit theorem affects eariier. This is due to the earlier diffusion of price in a 
discontinuous model, such as the jump model. In orderto obtain the same result as the Boltzmann model, it is necessary 
for the junnp model to use a distribution having a larger probability density in the region of lower eaming rate than the 

25 curve C1 0, thereby reducing the diffusion of prk^e. However, the resultant distribution of the daily eaming rate becomes 
quite different from the underlying assets. In this nrtanner, the jump model differs from the Boltzmann fundamentally. 
[0272] The Boltzmann model is not partbularly complicated, as compared with the jump model. This is tme from the 
comparison between the Boltzmann model and the simplest jump nrrodel, for instance, the Morton's complex jump 
model. The Morton's jump model uses a random number^ of the normal distribution, a random number of the Poisson 

30 distribution, a standard deviation a of the normal distributk>n. the average size k of jumps, and the probability X of 
occurrence of jumps per unit time. Thus, the Morton's jump model express stochastk; differential equation (45) of the 
underlying assets S using two probability density functions, that is, the Gaussian and the Poisson, and three param- 
eters. 

^ ^ = -Xkdt + a^7S + T| (46) 

[0273] In contrast, the Boltzmann nrtodel uses a single probability density function, that is, the Maxwell's distribution, 
and three parameters. The Boltzmann model is simpler using a less nunnber of probability density functions. 
^ [0274] The foregoing is the brief description of the related theory and the reality of option pricing. A dealing system 
for evaluating option prices based on the Boltzmann model will be explained hereafter. 

[0275] Rgs. 28 and 29 illustrate the structure of a dealing system 1 00 of the preferred embodiment of the present 
invention. The dealing system lOOconnprises an implied volatility cabulation unit 1 02, a Boltzmann model computation 
engine (BM M) 1 03, an implied volatility (IV) filter 1 04, and a dealing terminal 1 05. The implied volatility calculation unit 

^ 1 02 connmuncates with an external mark^ datatjase 101, and acquires nrtarket data to cakuilate impi led volatil ity. The 
Boltzmann nrtodel computation engine (BMM) 1 03 has a structure shown in Rg. 29, and conducts option price evaluation 
based on the Bottzmann model. The implied volatility (IV) filter 104 converts the option price provided from the BMM 
1 03 into the implied volatility (IV). The dealing terminal 1 05 functions as a graphical user interface (GUI), and displays 
necessary information. The dealing tenminat 1(^ also outputs data as, for example, hard copies, and inputs data to 

so the system 100. 

[0276] Rg. 29 illustrates the stmcture of the Boltzmann model connputation engine (BMM) 1 03. The BMM 1 03 com- 
prises an initial value setup unit 3, an evaluation condition setup unit 4, a Bottzmann model analysis unit 5, a graphk^ 
user interface (GUI) 105, a total cross^section^ochastic process setup unit 7, a velocity/direction distribution setup 
unit 8, and a random numt>er generator 9. The initial value setup unit 103 allows the initial values of at least a price, a 
^ price change rate, and a price change direction to supply to the BMM 1 03. The GU1 1 05 may be used in common with 
the GUI shown in Rg. 28. The BMM 1 03 is connected to an external market database 101 in order to take necessary 
market data in. 

[0277] The Boltzmann model analysis unit 5 includes an initializing unit 1 2, an initial value setting unit 1 3, a sampling 
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unrt 14, a price fluctuation simulator 15, a probability density calculation unit 16. rt also has a one-trial completion 
detector 17, an all-trial completion detector 18, a probability density editor 19, a price distribution calculating unit 20, 
and price converter 21 . The simulator 15 conducts calculations of price fluctuation based on the Boltzmann model. 
[0278] The invention will not be limited to the exemplified systm configuration, in which the system 1 00 is installed 

5 in a single computer in the physical nrteaning. For example, the system 1 00 can be a divided system, such as a client- 
server system for conducting operations separately. Preferably, the elements in the system correspond to programs 
for causing the associated elements to execute the operations, or perfonn the functions, indicated in the blocks in the 
drawings. Accordingly, the dealing syst^ 100 can be realized by installing a dealing program for causing a single 
computer with a communication function to perfonrn these functions. 

10 [Q279J The initial value setup unit 3 inputs Tq, Cq, gQ of equation (43) to the Boltznnann model analysis unit 5. These 
parameters are directed to the underlying assets of a stock price or a stock index to be evaluated, and are obtained 
from the past records. Preferably, the initial value setup unit 3 retrieves inf onmation on the stock price or the stock from 
the market database 101, and acquires the initial values of the price, the price change rate, and the price change 
direction from the retrieved information. The evaluation condition setup unit 4 supplies evaluation conditions to the 

IS Boltzmann model analysis unit 5. The evaluation conditions include the number of trials, a time band, a price range 
concerned, whch are required by the Boltzmann model analysis unit 5 to conduct essential calculations. 
[0280] The Boltzmann model analysis unit 5 is the center of the dealing system, and the most essential element. 
The stmcture of the Bottznnann model analysis unit 5 are almost the same as that of the price and risk evaluation 
system shown in Rg. 1 , but new elements, a price distribution cafeulation unit 20 and a price conversion unit 21 , are 

20 added. 

[0281] The price distribution calculation unit 20 cak:ulates the price distributk>n based on the price change probability 
density of the underlying assets edited by the probability density editor 19. 

[0282] The price converston unit 21 computes and outputs an option price, based on the price distribution calculated 
by the price distribution cak^jlation unit 20. 

25 [0283] The dealing terminal 105, whk:h functions as a GUI, outputs the intermediate operation results and the final 
results of the process. The deaUng terminal 105 also outputs the price distribution of the option to be evaluated. This 
terminal 1 05 has an input function for allowirtg the user to input data by a pointing device, such as a keyboard or a 
mouse. The dealing terminal 1 05 has an output furYCtions, such as displaying information on the monitorscreen, printing 
out as hard copi^, transferring data to other systems via a networic, and writing in to the memories, etc. 

30 10284] The market database 101 stores information as to option products. In this context, the tenrn "database" in- 
cludes the data systennatically stored in the database, means for retrieving the data, and the hardware storing and 
managing these. 

[0285] The market datat^ase 101 may be dedcated to this system, or altematively, it may be an existing extemal 
database, if there are any. 

35 [0286] An evaluation method of stock index option price carried out by the above-explained system will now be 
descrit>ed. 

[0287] Rg. 30 illustrates six steps A1 through AS. In the evaluation of stock index option price, many transaction 
data are generally utilized. The past records including such transaction data are accumulated in the market database 
101. In step A1, the total cross-section/stochastk: process setup unit 7 cateulates an implied volatility using the past 

40 records stored in the datat>ase 101 . 

[02^] The steps A2 and A3 illustrated in the left half of Rg. 30 show the conventional procedure for the sake of 
comparison. In A2, a general trend of a votatility matrix is detenmined from the implied volatility obtained in A1 , based 
on the experierrces and intuitk>n, or othenwise, on the simple path average or a recun^ent model. Then, in A3, a volatility 
matrix is determined. The conventional technique required ariliitrary judgments in step of A2. 

^ [0289] In contrast, thedeaUng system 1 00 of the prefen^ embodiment receives the implied volatility at the Boltzmann 
model analysis unit 5 in step A4, and delenmines the temperature parameters (i.e., the three coefficients Tq, Cq, in 
equation (43)) of the Boltzmann model. These parameters must be determined so that the temperature of the Boltzmann 
model coincides with the implied volatility. 

[0290] In step A4, it is determined wither the Boltzmann model temperature agrees with the implied volatility. If they 
50 agree with each other, then the process proceeds to st^ A5. in which it is detemiined whether the outcome of the 
Boltznnann model agrees with the daily earning rate of the underlying assets. If they do not agree with each other, the 
process returns to step A4, and the parameters are reselected. If they agree with each other, then the Boltzmann model 
is compared with the market-dependency of the underlying assets in step AS. 

[0291 ] If, in step of A6, the Boltzmann model agrees with the market-dq;>endency of the underiying assets, the proo- 
fs ess proceeds to step A3, in whch a volatility matrix is detemnined based on the outcome of the Boltzmann model. In 
reality, it is rare that an apparent market-dependency is observed. Accordingly, the flow from A5 to A6 becomes the 
final decision in many cases. If an apparent market-dependency is observed, and if the Boltzmann model tennperature 
does not agree with this nrtarket-dependency, then the process returns to step A4 in order to reselect the parameters. 
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[0292] The Boltzmann model can not always explain the real market very well. Even if a volatility matrix can be 
described very well in step A4, it may be contradict with the daily eaming rate. In such a case, the agreement between 
the Boltzmann model and the daily eaming rate may be given up, and the process proceeds to step A3. If the Boltzmann 
model does not agree with the implied volatility in step A4, it means that the actual market exceeds the limit of the 
5 Bolt^ann model. In this case, the process returns to the conventional step A2, and the evaluation is left to the judge- 
ment of dealers. 

[0293] Next, how the dealing system 1 00 evaluates the option price of an individual stock will be explained. 
[0294] It is observed that the daily eaming rates for the option estimated from the past records using the Boltzmann 
model agree with the reproduced option price based on the underlying assets through the Boltzmann model well. This 

10 fact means that the Boltznrtann model has capability for estimating the option prices based on the past records of the 
underiying assets, without past option-transaction records. Accordingly, the method of the present invention Is the most 
efficient method for evaluating the option price of an individual stock having little transaction record at present time. 
[02^ Basteally, the Boltzmann model uses the same technique for evaluating the option price of an individual asset 
as that for evaluating the option price of a stock index. In other words, the option prices of individual asset can also be 

IS evaluated by determining the three coeffictents Tq, Cq, go of equation (43) for each asset. 

[0296] Rg. 31 illustrates the applk»tk>n of the tendencies reproduced by equation (43) to individual asset. Rg. 31 
illustFEttes the temperatures T of various assets listed in the Tokyo Stock Exchange. Along the horizontal axis, vanous 
industries (constructions, foods, chemicals, steels, electric appliances, finances, Eind servk^es) are arranged in the 
increasing-order of the historical volatility of stock pnce. Solid lines C21 , C22, ... , C27 in Rg. 31 represent the temper- 

20 ature T converted from the historical volatility. In the drawing, the darkened circles (symbol #) denote the temperature 
of the eaming-rate distribution of the cunrent day with the daily eaming rate of within 5% of the previous day. The cross 
marks (x) denote the temperatures virrth the daily eaming rate of between 5% to 10% of the previous day. The white 
squares (□) denotes the temperature with the daily eaming rate of between 10% to 15% of the previous day. 
[0297] Fig. 31 reveals that the temperatures of three ranges are almost proportional to the corresponding historical 

25 volatility such as C21 , C22, C27. It is also found that the temperature T increases as the daily eaming rate of the 
previous day becomes larger. Conceming the distn'butk)ns of the three synntiols, the separatk>n between the group of 
the white squares and the group of the darkened circles (•) is larger than the separatbn between the group of the 
white squares (□) and the cross marks (x). This means that as the daily eaming rate of the previous day increases, 
not only the temperature, but also the tenrtperature rising rate increases. These observations suggest that temperature 

30 T has the quadratk: form with the prevk>us day's returns as expressed in equation (43). 

[0298] Figs. 32 and 33 illustrate examples of option price evaluation for individual stocks. Fig. 32 shows examples 
of evaluation for call optk)ns. The horizontal axis denotes the ratio of exercise price to the underlying price (i.e., the 
price of the underiying assets). The vertical axis denotes the ratio of the call option price to the underiying price. Fig. 
33 shows examples of evaluation for put options. The horizontal axis denotes the ratio of exercise price to the underiying 

35 price. The vertical axis denotes the ratio of the price of put option to the underiying price. In Figs. 32 and 33, the solid 
lines C31 , C41 represent the evaluated results of the Boltzmann model with the term to the maturity of 20 days. The 
one-dot broken lines and C42 represent the evaluated results of the Boltzmann model with the tenm to the maturity 
of 40 days. The dashed lines and C43, and the dotted lines C34 and C44 represent the results based on the Black- 
Sholes equation of the con^sponding terms. 

^ [0299] Here, temperature T is given in the foltowing equation (46) from Rg. 31 under the assumption that an under- 
tying asset concerned has historical volatility with the magnitude of about 70%. 



[0300] The magnitude at>out 70% of historical volatility is aIrTU>st doubled with the volatility of the stock index like 
option of Nikkei 225 stock average. Although this value is slightly larger than the historical volatility for many stocks, 
so the magnitude is still in the range of realistk: values. 

[0301 ] Rgs. 34 arKt 35 illustrate the innplied volatility for a call option and a put option , respectively. These two graphs 
exhibit smile curves and the term structures, as the smile curves and the term structure for the example of the implied 
volatility for the stock index shown in Fig. 27. 

[0302] Fig. 36 shows the operation flow in evaluating the option price for the individual stock executed by the dealing 
55 system 1 00. The left half of the drawing the conventional fkiw, as in Rg. 30. 

[0303] In general, the real records of the option transaction for yarwus stocks are insufficient to determine the his- 
torical volatility as the reliable infonrtation for option dealing. The option transactions with adequate records for various 
assets are handled with the analogous ways to the option dealing of a stock index, like the Nikkei 225 stock average. 
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In the following explanation, the flow of Rg. 36 is applied mainly to the option dealing with little records for inactive 
transactions. 



[0304] With the conventional technique, the smile and the term struc^re of the volatility are intended from the expe- 
riences in evaluation of option prices for stock Indexes, based on the historical volatility of the associated asset, in step 
B1 . Then, in step B2, a volatility matrix is detemnined based on the smile and the term structure. 
[0305] In contrast, the dealing system 100 of the prefen^ed enobodiment determines the temperature parameters of 
the Boltzmann model for this asset in step B3. Then, it is confirmed whether the daily earning rate of the Bottzmann 
model agrees with that of the underlying assets. If the daily earning rates do not agree with each other, the process 
returns to step B3 to reselect the paranr^ters. If they agree with each other within the spedf led accuracy, the process 
proceeds to step B5, in whic^ the agreement between the market-dependency of the Boltzmann model and the actual 
record is checked. When the both values agree with each other within the certain range, the process proceeds to step 
B2 to determine a volatility nrmtrix. When there is dtsagreenrYent in step B5, the process returns to step B3 to reselect 
the parameters, and to repeat the above-mentioned procedure. In reality, an apparent market-dependency is rarely 
observed, and therefore, the step B4 may be the final determination to proceed to step B2 (as indk:ated by the dashed 
arrow). 

[0306] Next, how the dealing system 100 evaluates an option price matching with the historical information, based 
on the Boltzmann model, will be explained. It has already been mentk)ned that the Boltzmann model Is capable of 
evaluating the smile and the temn structure of the implied volatility, while keeping consistency with the historical infor- 
mation of the undertying assets for reprodi^ng the derivatives. This feature is very advantageous, because the basis 
of the option price is dearly indcated when conducting transaction. In the current status for optbn pricing in the real 
market, rationale for pricing is not necessarily required. Espedally when the transaction is carried out within the full 
scope of setf-resportsibility, large losses due to mispricing from the dealer's judgement is attributed to a problem of the 
person concemed. On the other hand, In the consulting business for transaction or option pricing, it is required that 
pricing options must be determined with the reasonable bases, other than long experiences and intuition k>asis. 
[0307] Since the uncertainty is high in the market, arbitrary judgment relying upon experiences and intuitions will 
never completely disappear. However, if such arbitrary juc^ment is rationalized by other adequate information, the 
judgement becomes an action based on the rationale, and is not mere manual judgement any longer. The modem 
financial, engineering stands for foundation that pricing derivatives are uniquely detemrtined by the price of the under- 
lying asset. Therefore, the match between the Boltzmann model and the historical infomnation of the underlying assets 
can be a strong basis of rationality of price evaluation. 

[0308] It is recognized that prk^ng options by the jump model or volatility matrix model, previously described, cannot 
fully cover wider range of the historical informatk>n. In these cases, it is difficult to explain reasonable sources for the 
deviation from the reality in the option market. Since such nrtodels are not perfect, option pricing needs the new model 
whbh can reproduce option price from the historical informatk>n from the underlying assets. However, the Boltzmann 
hnodet is advantageous because it keeps consistericy with the historbal Inf omnation on the undertying assets in principle 
with a very little deviatk)n, and in addition, there are a relatively few occasions that require deviation. Even if such 
deviation is observed, dealers can discuss with the interested parties in advance. Therefore, a situation in whch mis- 
judgement by a particular dealer or a consultant amy lead to a huge amount of losses can be avoided. 
[0309] The dealing system 100 must evaluate risk parameters expressed by equations (47) through (51) as a risk 
hedge when evaluating a option prk^e. 



dC 
dS 



(47) 
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here C is the option price. S is the underlying price, r s the non-risky interest rate, t is the temri to the maturity, and a 
Is the volatility. It is known that if the underlying assets are traded in proportion to these risk parameters, the pnce 
fluctuatk)n of the underlying assets can be carrceled out in principle. 

[0310] These risk parameters are the differentials of the option price. The Boltzmann model uses the Monte Cario 
method for numerical analysis as a prerequisite; however, the Monte Cario nr^thod has a drawback of taking compu- 
tatk)n time for evaluating the differentials. For instance, when calculating the risk parameter e of a call option price C 
strictly using the Monte Cario method, an infinitesimal change 5t in temn xto the maturity is set, and equatk)n (52) is 
computed. 



[0312] The integral of the right-hand skie of equation (53) Is computed by the Monte Cario method. The change in 
T of the input variable is very small, and therefore, the difference t)etween the integral of the first tenm and the integral 
of the second term of the right-hand side of equation (53) is also small. 

[0313] Since the computation result of the Monte Cario method varies within the range of the statistics error, the 
amount of computation must be increased in order to reduce the statistic error. In general, the statistk: en^or is inversely 
proportional to the square of the computation amount Accordingly, if the change is very small, a significant difference 
can not be detected unless a vast of time is spent for computation. 

[0314] This problem applied to not only the financial Monte Cario, but also the general Monte Carlo, and an ultimate 
solution has not been found yet. The neutron transport Monte Cario method uses a perturt)atlon Monte Cario, whteh 
simulates only an infinitesimal change. However, even with this method, some approximation is required, and the 
advantage of the Monte Cario method, that is, the exactness may be damaged. 

[0315] Considering the fact that the Monte Cario method can not conduct a rigorous simulation of an infinitesimal 
change at present, It is not practical to stick to the exact solution. Under a circumstance in whk^ the option price can 
be explained by the implied volatility, separation b^een the Boltzmann model and the Biack-Sholes model is small. 
Because the first-order or the second-order differentials, such as risk paranrteters, do not depend on models very much, 
risk parameters can be evaluated at a practrcalty sufTictent accuracy by inputting the implied volatility calculated back 
from the option price of the Boltzmann model into the Black-Sholes risk parameter evaluation equation. 
[0316] To be more precise, the vo^'lity a of risk parameter evaluation equations (54) through (56) based on the 
Black-Sholes equatk>n is replaced by the implied volatility that coincides with the option price of the Boltzmann model. 
Then, the Black-Sholes risk parameters are used as those of the Boltzmann rTK>del. 




dr 

Call(K,T 6t) - Can(K,T) 
6r 



(52) 



[031 1] The numerator of equation (52) becomes equation (53) from equation (36). 



Call(K,T + 6x) - Can(K,T) = 

e-"( e-'* JdS(S- K)P(S,T + dT)-J'dS(S-K)P(S,T)) 



(53) 




(54) 
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here, Erf (x) is an equation d^ned by 



(58) 



40 



Erf(x)--^-exp 



(59) 



45 [0317] The dealing system 100 may have an altemalive function, that is, a function for generating a table from the 
probability density functions evaluated by the Boltzmann model, and for computing the option price based on the 
Riemann sum of vectors, instead of on the recalculation of the Monte Carlo method. 

[0318] Implied volatility is a volatility of the Black-Sholes equation computed backward so that the option price eval- 
uated by the Boltzmann nrtodei agrees with the option price evaluated by the Black-Sholes (BS) equation. In the Boltz- 
so mann model, the option pnce is expressed by equations (36) and (37), and the integrals of these equations are obtained 
using the Monte Carlo method. 

[0319] If the probabilSy density function P (S, t) does not change largely with respect to S, that is, if a general pnce 
di^bution can be applied, then it is effective and useful to make a table of probability density functions with r^pect 
to S, and to make approximation by the series of equations 60 and 61 . In this case, the approximated result becomes 
55 quite dose to the rigorous evaluation result of the Monte Carfo method. 
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Can{K)»e-'' 2s!rK^Si P(Si,T) (Si- K) (60) 

5 

Put(K)-e-" Ssil^ASi P(Si.T) (K-Si) (61) 

10 

here the probability density function P (Si, t) Is expressed by equation (62), and is evaluated by the Monte Carlo method. 

f'"M^dS P(S,T) 

P(Si,T). '"\si <«^> 

20 [0320] By storing the probability density functions of equation (62) as a table with a practically sufHctent small Asi, 
equations (60) and (61 ) can be computed at a high speed because these two equations are merely the sum of product 
(or the Riemann sum) of vectors. 

[0321] Next, the operation process by the dealing system 100 will be explained. Rg. 37 illustrates a sub-screen (or 
a window) of the dealing temninal 105, which displays the detailed track of stock index in a continuous session. Fig. 
25 38 illustrates a sub-screen of the dealing terminal 1 05, whk^h displays a table of implied volatility and market prrces of 
each delivery month and each exercise price of stock index option together with the stock index as the underlying asset 
prices. 

[0322] Rgs. 39 (a) and (b) also illustrate sub^screens of the dealing terminal 105, which display the information 
contained in the table as graphs. Rg. 39(a) is a graph of implied volatility as a function of exercise price, which exhibits 

30 a so-called smile curve. Rg. 39(b) is a graph of option price of each delivery month as a function of exercise price. 
[0323] This ctealing system 1 00 carries out the operation flow shown in Rg. 40. In the nonnat maricet state, the system 
takes the market data 101 in (step SOS), and conducts rough calculation (step S10). The rough calculation result is 
displayed on the sub-screen shown in Rg. 37 (step S15), thereby improving the display speed. 
[0324] If £in abrupt change occurs to the underiying assets as indicated by the dashed circle "a" in Rg. 37, the user 

35 can switch the sub-screen (or the window) shown in Fig. 41(a) to a sub-screen of detailed information shown in Fig. 
41 (b). To be more precise, the user enlarges the window by dragging a desired area 111 (corresponding to K4-K5) 
shown in Rg. 41(a) in a diagonal direction using a mouse or other pointing devk» in a diagonal direction (steps S 20 
and S25). Then, additional columns are produced corresponding to the dragging amount, as shown in Fig. 42(b). The 
additional input data, that is, virtual exercise prices K451 , k452, artd K453 are autonnatically transferred to the BMM 

40 1 03 for interpolation (steps S30-1 and S30-2). 

[0325] The interpolation result, that is, the implied volatility and the option exercise prices corresponding to the abrupt 
change "a" are retumed from the BMM 1 13 to the dealing terminal 1 05 (step S30-3). Then, the scale of the sub-screen 
of F^, 38 and 42(a) and the scale of the sub-screen of Rgs. 39 and 41 (a) are refreshed (step S 35-1 ) to display the 
detailed information. The interpolated results (I.e., the ctetaited information) are connected into a line, and displayed 

45 on the sub-screen (step S35-2). F^. 41 (b) shows the interpolated detailed implied volatility C65, in comparison wrtii 
the rough calculation result C64. 

[0326] Next, the computation carried out by the BMM 1 03, whk:h is a theoretical computation server, will be explained 
with reference to Rg. 43. 

[0327] As has already been mentioned, implied volatility (IV) is a volatility of the Black-Sholes (BS) equation calcu- 
so lated backward in such a manner that the option price evaluated by the Boltzmann model agrees with the option price 
provided by the BS equation. 

[0328] In the Boltzmann model, the option prices of call option and put option are described by equations (36) and 
(37), and the integrals of these equations are obtained using the Monte Cario method. If the probability density function 
P (S. x) in equatk>ns (36) and (37) does not exhit>it an extremely big change, that is, if the ordinary price distribution is 
S5 applicable to the evaluation, then, the prot>ability density function when it is not largely varied, i.e., when the normal 
price distribution can be applied, then the probability density functions P (S.-c) are arranged in a table with respect to 
S. By making approximation using the series of equations (60) and (61), a consklerably accurate result dose to the 
strict evaluation result by the Monte Cario method can be obtained. 
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[0329] Here, the probability density function P (Si, x) contained in equations (60) and (61 ) is expressed by equation 
(62), which is evaluated by the Monte Carlo method. By taking AS\ sufTtctentty and practically small, and by storing the 
prot>abtlity density function of equation (62) as a table, equations (60) and (61) become merely Ri^ann sums of 
vectors, and can be computed at a high speed. 

5 [0330] Rg. 43 schematically illustrates the series of equation (60). A smooth curve C61 in Fig. 43 is the genuine 
probability density P (S, t). The histogram C62 is the tabled probability density P (Si, t). The Monte Carlo integral of 
the product of the genuine probability density P (S, t) and the linear line C63 with a slope 1 and starting from the point 
of S K is the strict solution of equation (36). If the user requests an interruption for evaluating a detailed result (YES 
in step S20 of Rg. 40), the rigorous solution of equation (36) is computed. On the other hand, the integral of the product 

10 of the linear line C63 and the histogram C62 is an approximation expressed by equation (60). In the normal market 
state (NO in step S20), approximation is computed as a rough evaluation. 

[0331] In this manner, if detailed evaluatk)n is desired, stricter result is obtained and displayed. In the graph of Fig. 
41 (b), a symbol C indicates the separation between the ordinary rough cak^lation 064 and the detailed evaluation 065. 
[0332] With this arrangement, the evaluation result is displayed pronrptly based on the rough cakxilation in the or- 
dinary state. Under a user request, the window or the sub-screen is enlarged to display the detailed evaluation, thereby 
allowing the user to visually know the market change qukHdy. 

[0333] tMext, a technique for developing a structured bond or an exotk: option by the dealing system 100 shown in 
Rg. 28 will be explained. Such a development can be realized by obtEiining and displaying volatilities of arbitrary multiple 
terms to evaluate the term structure of the volatifity independent of the market Fig. 44 illustrates the operation flow of 
20 evaluatk>n of multi-term volatility. The same steps as those in the flowchart of Rg. 40 are denoted by the same numerical 
references. 

p>334] In the ordinary state, the dealing terminal 105 displays the market on the sub-screen based on the rough 
calculation result, as shown in Rgs. 37 through 39. If the user wishes to evaluate volatilities of multiple terms, a tenti 
setting window 200 shown in Fig. 45 is called. In order to select the evaluation period for a selected option, a starting 

25 date 201 , the maturity 202, and the evaluation interval^ are input 

[0335] The input information does not have to be those values set in the market For example, the options in the first 
and second delivery months circled with 'a" in Rg. 46, whk:h shows the implied volatility and the market prices, are 
traded in the market However, there is no transaction of the option of the m*^ delivery month, which is circled with "b". 
[0336] If the user requests an interruption (step 820* in Rg. 44) for evaluation of the m*^ delivery month, the user 

30 inputs the final date of the m^ delivery month in the item 202 on the screen 200 with the selection of monthly evaluation 
(step S25'). Then, the inputted infomrtation is automatically transfen-ed to the BMM 103 for interpolation (step S30'; 
steps S30-V and S30-2^ upon hitting "enter". 

[0337] The BMM 1 03 carries out the connputation based on Boltzmann model. The dealing tenninal 1 05 receives the 
computation result, more specifically, the implied volatility and option price "b" of the m^ delivery month (Step S30^'). 
35 Then, the scale of the sub-screen shown in Figs. 38 and 39 are refreshed, and the interpolation results are connected 
into a line and displayed on the refreshed screen (steps S35-1 , S35-2; 835). 

[0338] Rg. 47 illustrates an example of the interpolation result, in whk^h the curves 071 and 072 represent the 
volatility of the optbn existing In the market, and a curve 073 represents the volatility of the option that does not exist 
in the market. 

^ [0339] In this manner, obtaining and displaying the volatility of arbitrary multiple temns allows the tenn structure of 
the optk>n volatility that does not exist In the market to be evaluated. This arrangement can improve the efficiency of 
developing a structured bond or an exotic option. 

[0340] Next, a technique for displaying the behavior of the term structure of the implied volatility (IV) at the money 
(AT^) with fading animations will be explained with reference to Rgs. 48 through 52. The fading animation is also 
^ carried out by the dealing system 100. 

p)341] The flowchart in Rg. 48 Includes the facfing steps, whch are basbally added to the operation flow shown in 
Rg. 40. The same steps as those in Rg. 48 are denoted by the same numerical references. The newly added function 
is the ^ding display in step S150. 

[0342] In the ordinary state, the market is displayed on the sub-screen of the dealing temninal 1 05, as shown in Figs 
50 37 through 39 based on the rough cakujIaUon (step S10). If an abrupt change occurs in the underiying assets (as 
Indbated by "a" in Fig. 37), the user can select a desired area 110 in the sub-screen of Fig. 39 to took into the state. 
The user can enlarge the selected area for detailed evaluation, as shown in Rg. 41 (a), using a mouse (YES in step S20). 
[0343] In response to this user intenruption, the necessary data provided by the user in the enlarged screen is auto- 
matically trartsferred to the BMM engine 1 03. The simulation result by the BMM engine 1 03 is retumed to the dealing 
55 terminal 1 05, and displayed on the sub-screen, as shown in Rg. 50. The column KR shown in Fig. 50 represents virtual 
realtime ATM interpolated between the discrete values of the exercise price bands set in the maricet. The virtual ATM 
is assumed to be equal to the real-time underlying price. 

[0344] Since the risk index and the option price change most signifk^antty near ATM, it.is very important for dealers 
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and traders to observe the market change and the term structure in this area. 

[0345] With the dealing systemlOO of the present invention, even if the underlying price nnoves between discrete 
exercise price bands set in the market, virtual ATM is interpolated as indicated by KR in Rg. 50. Dealers can flexibly 
evaluate the term stmcture of the implied volatility based on the interpolation. 

5 [0346] However, a problem is that the anrraunt of inf omnation displayed on the graph of Rg. 39(b) inevitably increases. 
For example, if the user wishes to evaluate up to the sixth delivery month, six option prices must be displayed on the 
graph. In this case, it is desirable for the user to know these evaluation results conrectty and promptly at a glance. 
[0347] To respond to this demand, the behavior of the tenm structure of the implied volatility is displayed in fading 
animations, as shown in Fig. 52, based on the operation flow shown in Rg. 49. In Fig. 52(a), the solid line draws only 

10 the implied volatility of the a1 ^ delivery month (step SI 504). After a time lapse of a weight level w, the implied volatility 
of the a1^ delivery month is displayed by the dashed line, and the implied volatility of the a2'^ delivery month is 
displayed tyy the solid line, as shown in Rg. 52(b) (steps SI 50-5, S1 50-6). This step is repeated until the ae^ delivery 
month, as shown in Rg. 52(d) (steps S150-7. S150-8, S150-e, S150-(et-1)). 

[0348] The weight level w can be adjusted by typing the plus (+) key or the minus (-) key of the keyboard, as shown 
IS in steps SI 50-1 through S150-3 in Rg. 49. In response to the key manipulation, an intenruption routine is activated, 
and this interruption is reflected to the fading display shown in Fig. 52 at once. The interruption (i.e., adjusting the 
weight level) may be carried out by using a mouse or a pointing device, instead of keyboarding. The computation result 
of Rg. 50, which was retumed from the BMM engine 103, is display as a grEiph, as shown in Fig. 51 . In Rg. 51 , KR 
denotes the virtual ATM, the white triangles denote the output data from the BMM engine 1 03, and the cross marks 
20 dertote the market cteta. 

[0349] Displaying the behavior of the term structure of the implied volatility (IV) in ATM by fading animation can 
protect the user from mispricing due to too sensitive reaction to the market change. 

[0350] Fig. 53 is a flowchart including a fading function, whk^ is basically added to the operation flow of Fig. 44. The 
same steps as those in the flowchart of Rg. 44 are denoted by the same numerical references. The newly added step 

25 is the display processing of step SI 50*. 

[0351 ] In the flowchart of Rg. 53, the market color associated with market activity is displayed on the sut>-screen of 
the dealing tenminal 105 in the ordinary state based on the rough caknilation result, as shown in Figs. 37 through 39 
(step S10). If the user wishes to evaluate the votatilrty of arbitrary multiple temns simultaneously, the user inputs the 
rtenr^ of staring date 201, the maturity 202, and the evaluation inter/al 203 through the temi setting page shown in 

30 Rg. 45 (step S20') to set the tenm of a desired option. 

[0352] In response to this user interruption, the data inputted by the user is automatically transferred to the BMM 
engine 103. and the computation result shown in Rg. 46 is retumed (steps S25' and S30'). This computation result is 
displayed as a graph, as shown in Rg. 47 (step S35). 

[0353] If the user designates the fading mode, mutti^emn volatility is displayed on the screen by fading animation 

35 based on the operation ftow shown in Fig. 49. 

[0354] By displaying the animated behavior of the tenm structure of the implied volatility of arbitrary multiple terms 
in a fading manner, undesirable mispricing due to excessive reac^'on to the market change can be prevented. 
[0355] Another significant function of the dealing system 100 shown in Rg. 28 is to allow dealers to set positions 
and to automatically order desired financial products. This function will be explained with reference to Figs. 54 through 

40 56. In Rg. 54, the same steps as those in the flow of Rg. 40 are denoted by the same numerical references. 

[0356] In the ordinary state, the sut>-screen of the dealing tenminal displays the market color associated with market 
activity based on the rough-computation as shown in Figs. 37 through 39 (steps S05 through SI 5), and the user can 
visually refer to the appropriate standard of the advanced model. 

[0357] When the user d^res to order transactions, the user requests an intenruptbn at step S40, and drags the 
45 areas L, M and N shown in Rg. 55 In turn using a mouse (steps S40 and S45). If the user first designates the area L, 
thetransaction input page 210 shown in Rg. 56 is started. The page 210 exhibits the attributes of the designated 
position, such as exercise prices, call/put option, implied volatility, and limit price. The user can select in this page 
purchase or sell, together with the desired numt)er of bOts. 'Upon hitting the 'enter" key, the designated order is stored 
in the order database 130 (step S50). The other areas M and N are handled in the same manner. 
50 [0358] After the order was made, it is. checked every t seconds in real time whether the market has regressed to the 
target area (steps S60, SI 0, SI 5, 820 and S40). 

[0359] If the market has regressed to the target area (YES in step 860 of Fig. 54), the user's order is transferred to 
the market immediately and automatically (steps S65 and S70). If the position in the target area (e.g., area N) has 
been contracted as a trade date (step S76), then area N is eliminated from the page (step S75). 
55 [0360] This operation flow allows the user (or the dealer) to visually check the appropriate standard of the advanced 
model, and to nrtake orders timely in an autonnatk: manner. 

[0361] . This automatic ordering function based on a selected position may be combined with a fading function. In this 
case, the behavior of the term structure of the implied volatility ATM is displayed as a fading animation on the dealing 
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terminaJ 105. Figs. 57 through 61 illustrate the operation flow of the fading function (wnbined with the autom 
function. The same steps as those in other flowcharts are denoted by the same numerical references. 
[0362] In the ordinary time, the market color is displayed on the sub-screen of the dealing terminal 1 05 based on the 
rough computation result, as shown in Rgs. 37 through 39 (steps 505, and S10). 

5 [0363] Since the risk inctex and the option prk:e change significantly near ATM, it is very important for dealers and 
traders to observe the market change and the term structure. With the dealing system of the present invention, even 
If the underlying prrce changes between discrete exercise price bands set in the mari(et on the sub-screen shown in 
Fig. 37. ATM is interpolated as indicated as symbol KR in Fig. 62. The term structure of the implied volatility is then 
displayed by fading animatton (steps S150-1 through S150-(e+1)). 

10 [0364] Based on the infonnation displayed on the dealing terminal, the user dicks a desired position (for example, 
area L) in the displayed page of Rg. 62 (steps S40, S45, S50). It is automatically checked every t seconds whether 
the market has regressed to the target area in step S60. If the market has regressed to the target area, transaction 
order is immediately transfen^ed the market in an automata manner (steps S65, S70). 

[0365] In this manner, the dealing system 1 00 of the preferred embodiment allows the dealer to visually check the 
IS appropriate standard of the advanced model with a fading animation, and to select a desired position to make a trans- 
action order timely and automatically. Because the behavior of the term structure of the implied volatility (IV) in ATM 
is displayed by animation in a fading manner, the dealer can avoid misprk;ing due to an excessive reaction to the 
market change. 

[0366] Next, an automatic warning function of the dealing system 100 will be explained with reference to Figs. 63 
20 and 64. In Fig. 63, the same steps as those in the ftowcharts of Rgs. 57 through 61 are denoted by the same numerical 
references. 

[0367] In the ordinary state, the market is displayed on the sub-screen of the dealing terminal 1 05 based on rough 
computation, as shown in Rgs. 37 through 39 (steps S05, S10, S150), and the user can visually refer the appropriate 
standard of the advanced model. 
25 [0368] In the sub-screen shown in Fig. 64, curve CI 01 represents the output of the Bottzmann model engine 103, 
curve CI 02 represents the predicted regression of the market, and curve C 1 03 represents the implied volatility of the 
market. If the user wishes to make a transaction order at position "a" in Rg. 64, the user dk^ks the position to designate 
this area for a transaction order in step S80. 

[0369] At this point of time, the user can select the automatic warning function for automatically giving a caution for 
30 the selected position "a" in response to the market change. To be more precise, the user requests an interruption in 
step S80, and selects a warning area "b" on the sub-screen shown in Fig. 64 in step S85. The selected waming area 
corresponds to the market area for which the user needs waming. the designated waming information is stored in the 
waming area database 131. 

[0370] If the user selects multiple positions, and wishes to set a warning area for each position individually, the user 

35 repeats the intermption as many time as the user wishes. 

[0371] As shown in Rg. 64, the selected ordering position "a" corresponds to the option of the exercise price K2, 
and the waming area for this poston covers both the exercise prices K2 and K3. This is one of the significant features 
of this waming function. The user can flexibly take exercise prices (i.e., option issues) of other than the designated 
position as a waming corrdition. 

40 [0372] After the waming area has been selected, it is automatk^ly checked in step S90 whether the market has 
entered the waming area V. This determination is conducted by checking both the waming area database 131 and 
the market datakjase 1 01 in real tirrie every t seconds. If the prediction of the maricet fails and the real market represented 
by curve CI 04 enters the waming area 'b", then the position "a" at which a transaction order is requested is immediately 
fibkered, as indicated by the hatched area *d" in Rg. 64 to caution the user (step S95). 

^ [0373] This waming function enables a risk manager to carry out appropriate risk management because a caution 
is auton^tically issued wh^ ttie nrmrk^ has entered the waming area designated by the dealer. An extra function for 
suggesting an altemative position may be ac^ied to the above-described waming function. In this case, the operation 
flow shown in Fig. 65 is inserted aft^ step S85 of Rg. 63, as indicated by symbol "D". If the predk^k>n of the existing 
ordering position fails, and if the real maricet enters the waming area "b" as indicated by curve CI 04, then an altemative 

50 position area "e" is suggested, as shown in Rg. 66, in addition to giving a caution. 

[0374] Prior to suggesting the atterrtatlve position area, the sub^screen of the dealing system has been displaying a 
rough computation result, as shown in Rgs. 37 through 39, and an ordering position "a" and the conresponding waming 
area "b" are designated on the sub-screen, as shown in Rg. 66. 

[0375] Now, if a change occurs to the real market, and if the market index CI 04 of the real market enters the waming 
55 area "b", the position "a" is inmmediately flickered to give a caution to the user. Then, in step S86 shown in Rg. 65, an 
altemative position area 'e* is cak^jlated based on the ordering data 131 and the maricet data 101. The altemative 
position area fails into an opposite positk>n for cx^mpensating the gap between the predk^ed position and the real 
market represented by CI 04. 
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[0376] The alternative position "e" is immediately displayed on the dealing temnlnal 105 in step S88, as shown in 
Rg. 66. 

[0377] The user can again make a transaction order at a desired position T inside the alternative position area "e". 
[0378] By adding the function for suggesting an attemative position to the warning function, a vast amount of loss 
5 due to an excessive reaction to the market change can be effectivety avoided in advance, because an alternative 
position is automatically extracted together with a caution. 

[0379] As has been described above, the dealing system of the present invention has many significant advantages 
over the conventional technique based on a general theory of financial engineering, the application of which is very 
limited in a non-active option market. The dealing system of the present Invention has a Boltzmann model computation 

10 engine for applying a theory of nuclear reactor to the finandal field, which is capable of providing a theoretical price or 
a risk ind^ meaningful to dealers and traders through an interactive graphical interface in a flexible manner. Accord- 
ingly, the dealing system of the invention can deal with a big price change in the underlying assets in a flexible manner 
[0380] Although the present invention has been described based on examples of stock index option as the prefen-ed 
embodiments, the invention is not limited to these examples. The present invention is applicable to any option products, 

IS for example, indivkiual stock options and currency options, the underlying assets of which exhibit a behavior of the 
geometric Brownian motions. 

Claims 

20 

1 . A price and risk evaluation system for evaluating a price distribution or a risk distribution for a financial product or 
its derivatives, conrtprising: 

an initi'al value setter configured to input at least one of initial values of a price, a price change rate, and a 
price change direction for a financial product or its derivatives; 

an evaluation condition setter configured to input evaluation conditions including at least time steps and the 
number of trials; 

a Boltzmann model analyzer configured to receive at least one of the initial values from the initial value setter, 
receiving the evaluatk>n conditions from the evaluation condition setter, and repeating simulation of price fluc- 
tuation based on the Boltzmann model using a Monte Cario method within the range of the evaluation condition 
to obtain a price distrbution or a risk distribution for the financial product or its derivatives; 
a velocity/direction distribution setter configured to input probability distributions of the price, the price change 
rate, and the price change direction for the financial product or its derivatives into the Boltzmann model ana- 
lyzer; 

a random number generator for providing a series of random numbers used in the Boltzmann model analyzer; 
and 

a unit configured to output analysts results of the Boltzmann model analyzer. 

2. The price and risk evaluation system according to daim 1 , wherein: 

the initial value setter acquires the initial values of the price, the price change rate, and the price change 
direction for the finartdal product and the derivatives from a market database storing information on financial 
products or their cterivatives, and supplies them to the Boltzmann model analyzer means; and 
the velocity/direction distribution s^er receives past records for a financial product or its derivatives from the 
mart(et database, and then generates a probability density func^on with variables concerning the price, the 
prk^e change rate, the price change direction, and time steps for supplying the probability density function to 
the Boltzmann model analyzer. 

3. The price and risk evaluation system according to daim 1 , further comprising a total cross-section/stochastic proc- 
ess setter configured to supply information conceming setting a sampling time width in the simulation of price 
fluctuation to the Boltzmann mo<^ analyzer, wherein: 

the total cross-section/stochastic process setter acquires a price fluctuation frequency and a price change rate 
for the financial product or its derivatives from a market database storing information on financial products or 
S5 derivatives, and supplies ratios of the price fluctuation frequency to the price change rate into a total cross- 

section for a Boltznnann*s equation. 

4. The price and risk evaluation system according to daim 1 , wherein: 
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the velocity/direction distribution setter acquires past records for a financial product or its derivatives from a 
nnaricet database storing information on finandal products or derivatives, and estimates a distribution of the 
price change rate for the financial product or its derivatives using a Sigmoid function and its approximation 
form, and supplies the estimated distribution of the price change rate to the Boltzmann model analyzer. 

5 

5. The price and risic evaluation system according to daim 1 , wherein: 

the velocity/direction distribution setter acquires past records for a financial product or its derivatives from a 
maricet database storing infonnation on financial products or its derivatives, detenrtines a set of Sigmoid func- 
10 tion parameters with the price change rates using the past maricet data for a distribution of the price change 

rate, and supplies the price distribution to the Boltzmann model analyzer 

6. The price and risic evaluation system according to claim 1 , wherein: 

IS the velodty/direc:tion distribution setter acquires past records for a financial product or its cterivatives from a 

maricet database storing infonmation on financial products or their derivatives, estimates the probability distri- 
bution of the price c^hange direction for thefinandal product or its derivatives from the past records, and supplies 
the probability distribution to the Boltzmann mociel analyzer. 

20 7. The price and risk evaluation system according to c^aim 6, wherein: 

the velocity/direction distribution setter estimates the probabilrty distribution of the price change direction for 
the financial product or its derivatives taking into accx)unt a correlation between a probabilrty of price-up and 
a probability of price-down. 

25 

8. The price and risk evaluation system according to c^aim 1 , wherein: 

the velocity/direction distribution setter rec^eives past records of a financial product or its derivatives from a 
rrtaricet database storing infonrtation on financial products or their derivatives, and supplies the probability 
30 distn'bution to the Boltzmann model analyzerth rough generating the probability distributions taking into account 

a correlatk>n between a distribution of the price change rate and a distribution of the price change direction 
for the financial product or its derivatives. 

9. The price and risk evaluation system according to c:laim 1 . wherein: 

35 

the velo(^ty/directk>n distribution setter generates homogeneous probability distributions independent of the 
price, or heterogeneous probability distnl>utk)ns depenctent on the price, with regard to the probability distri- 
butions of a price change rate and a price change ciirection, and supplies the probability distributions to the 
Boltzmann nrtodel analyzer. 

40 

10. The price and risk evaluation system according to davn 1 , wherein: 

the Boltzmann mociel analyzer uses a linear Boltzmann model or a non-linear Boltzmann model in order to 
price a financial product and its derivatives, the linear Boltzmann mcxlel using a cross-section independent of 
^ prot>ability density or ftux for the financial prociuct or its derivatives in a Boltzmann equation, while the non- 

linear Boltzmann model using a c:rDss-sectk>n dependerrt on the probability density or the flux for the finandal 
product or its derivatives in the Boltzmann equation. 

11. The price and risk evaluation system ac(x>rcling to daim 1 , wherein: 

50 

the Boltzmann model analyzer provicies a pric^e distribution or a risk distribution for the finandal product or its 
derivatives by using the ftux defined in a Boltzmann equation as a bi-products of a probability density f uneven 
and a price c:hange rate per unit time for the financial product or its derivatives. 

^ 12. The price and risk evaluation system according to daim 1 , wherein: 

the Boltzmann model analyzer applies a track-length estimator nnethod using fluxes assoc:iated with the finan- 
c:ial products or their cierivatives, as a varianc^e redu:tion method for a Monte CaHo calculation, to evaluate a 
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probabilrty densrty at an arbitrarily specified point of tinne. 

13. The price and risk evaluation system according to daim 1 , wherein: 

the Boftznnann model analyzer makes use of a point detector method, whbh is effective in a neutron transport 
Monte Carlo simulation, as a variartce reduction as the variance reduction method to evaluate a price or a risk 
distribution at an infinitesimal price band or an infinitesimal time zone f orthe financial product or the derivatives, 
using all of or a part of the price fluctuation events observed in the simulation. 

14. The price and risk evaluation system according to daim 1 . wherein: 

the Boltzmann model analyzer calculates an adjoint probability density or an adjoint flux reduced from an 
adjoint Boltzmann equation for a price fluctuatbn for the f inandal product or its derivatives, and weights sam- 
pling values proportional to the adjoint probability density or the adjoint flux, thereby reducing variance. 

15. The price an6 risk evaluation system according to daim 1 , wherein: 

the velocity/direction distribution setter supplies a correlated probability density distribution to the Boltznrtann 
model analyzer, correlated probability densities being estimated taking conrelations between price change rate 
distributions and price change cfirection distrftyutions into account for any financial products and their deriva- 
tives. 

16. The price and risk evaluation system according to daim 1 , wherein: 

the Boltzmann analyzer conducts procedures for evaluating a price or a risk distribution for a finandal product, 
and for applying the Ito's theorem to obtain a price or a risk distribution for its derivatives. 

17. The price and risk evaluation system according to daim 1 , wherein: 

the Boltzmann analyzer indudes algorfthms conceming multiple processing in simulations and reducing the 
analysis results of the multiple processing to obtain the probability density distribution. 

18. A recording medium storing a program for controlling a computer, the program comprising the steps of: 

causing an initial value setter of the ctvnputer to input at least one initial values of a price, a price change rate, 
and a price change direction for a finandal product or its derivatives; 

causing an evaluation condition setter of the computer to input evaluatton conditions induding at least an 
evaluation time and the numt^ of trials; 

causing a Boltznrtann model analyzer of the conrtputer to receive saki at least one of the initial values from the 
initial value setter and the evaluation conditions from the evaluation condition setter, and to repeat simulations 
of price fluctuation using a Monte Cario method based on the Boltzmann model within the ranges of the eval- 
uation conditions to obtain a price distributbn or a risk dtetn'bution for the finandal product or its derivative; 
causing a velocity/direction distrft)ution setter to supply probability distributions of the price, the price change 
rate, and the price change direction for the finajK:tal product orthe derivatives; 

causing a random number generator of the computer to provide random numbers that are used for analysis 
in the Boltzmann model analyzer; and 

causing an output unit to provide artalysis results obtained by the Boltzmann analyzer. 

19. A dealing system comprising: 

an implied volatility connputation engine configured to evaluate an implied volatility based on market data; 
a Boltzmann model computation engine confrgured to evaluate an option price for a selected option produd 
based on the Boltzmann nnodel using the maricet data; 

a filter configured to convert the option price obtained from the Boltzmann model computation engine into a 
volatility of the Black-Sholes equation; and 

a dealing terminal configured to display the volatilrty of the Black-Sholes equation in comparison with the 
implied volatility cakxilated from the market data, or to display the option price caknilated by the Boltzmann 
model computatk)n engine in comparison with an option price in maricet. 
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20. The dealing system according to dalm 19, wherein the selected option product is a stock price index option. 

21 . The dealing system according to daim 1 9. wherein the selected option produd is an individual stock option. 

5 22. The dealing system according to daim 19, wherein the Boltzmann model computation engine has a unit configured 
to calculate an option price consistent to historical information from the market data. 

23. The dealing system according to daim 1 9, wherein the Boltzmann model computation engine has a converter that 
converts an option price, whch was obtained in a discrete manner with respect to an exercise price based on the 

10 Boltzmann model, into a volatility obtained from the Black-Sholes equation in order to obtain the option price and 

the risk parameter through interpolation with the Black-Sholes equation. 

24. The dealing system according to daim 1 9, wherein the Boltzmann model computation engine has a table generator 
that generates a table of a probability density function evaluated by the Boltzmann model, and caknjiates an option 

IS price using a sum of inner products of vedors. 

25. A dealing system comprising: 

a dealing terminal having a graphics user interface; 

20 a theoretical price and parameter computation engine configured to switch between a rough computation and 

a detailed computation for a theoretical option price and a parameter, the rough computation providing a the- 
oretical option price and a parameter for each exerdse price and for each delivery month set in the market in 
a normal maricet state, and the detailed computation providing detailed infomnation including theoretical prices 
and parameters for exercise prk^es and assodated delivery months that are not set in the market, but specified 

25 by a user; 

a theoretcal price interpolatbn unrT, and 

an interface configured to receive market data, 

the dealing system nomr^ly causing the dealing temninal to display a market activity based on the rough calculation, 
30 and to display a detailed computation result in a assodated price band when designated by the user. 

26. A dealing system comprising: 

a dealing terminal having a graphics user interface; 
35 a rough computation engine configured to compute a theoretical price and a parameter for each exerdse value 

and for each deUvery month set in the market; 

acomputation engine configured to compute theoretical option prices and parameters for arbitrary terms based 
on a Boltzmann model; 
a theoretical price interpolation unit; and 
40 an interface for receiving market data, 

the dealing system normally causing the dealing terminal to display the maricet adivity based on the rough com- 
putation result, and causes the dealing tenminal to display volatility parameters cak:ulated for the terms designated 
by a user. 

45 

27. A dealing system conrtprising: 

a dealing terminal having a graphics user interface; 

a rough computation engine configured to compute a theoretical price and a parameter for each exerdse price 
so and for each delivery month set in the market; 

a detailed computation engine configured to conpute detailed information induding theoretk:al prices and 

parameters for exerdse prices artd delivery months that are not set in the market; 

a theoretical price interpolation unit; 

a position setter; 
S5 an automatic orcter trarYsaction unit; and 

an interface configured to receive market data. 

the dealing system sending an automatic order signal when a stock index option price or an individual stock option 
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price reaches a specified automatic ordering price band. 

28. The dealing system according to claim 26 or 27, further comprising a fading processor, wherein: 

the fading processor causes the dealing terminal to display a term structure for a volatility that has been 
converted from an option price at the money obtained by the Boltzmann model with animation in a fading mode. 

2B. The dealing system according to claim 26 or 27, further comprising a risk limit setter, wherein: 

the risk limit setter alerts the user when the market enters a waming area. 

30. The dealing system according to claim 28, further comprising a risk limit setter, wherein: 

the risk limit setter alerts the user when the market enters a waming area. 

31. The dealing system according to claim 29 or 30, further comprising an alternative position selector, wherein: 

the alternative positk>n selector selects an alternative position, and causes the dealing tenminal to display the 
alternative position when the risk limit setter alerts the user. 

32. A conrtputer-readable recording medium storing a dealing program for use in a dealing system, the program com- 
prising the steps of: 

causing the dealing system to compute an Spiled volatility using market data supplied to the dealing system; 
causing the dealing system to compute an option price based on the Boltzmann model using the market data 
for a selected derivatives; 

causing the dealing system to convert the option price obtained from the Boltzmann model into a volatility of 
a Black-Sholes equation; and 

causing the dealing system to display the volatility of the Btack-Sholes equation in comparison with the implied 
volatility cak^lated from the market data, or to display the option price based on the Boltzmann model in 
comparison with an option price in a market 
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